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'^"^pThe  aim  of  this  study  is  to  provide  the  first  detailed 
look  at  the  morphology  of  the  ionospheric  F-region  near  a geo- 
magnetic latitude  of  60  degrees.  This  region  is  of  consider- 
able interest  because  it  lies  at  the  base  of  the  geomagnetic 
field  line  which  forms  the  average  boundary  of  the  plasmasphere 
the  geomagnetic  field-aligned  envelope  of  co-rotating  plasma  • — 
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which  extends  to  approximately  four  earth  radii  in  the  equator- 
ial plane  (L  = 4) . 

The  observational  parameters  used  in  this  study  monitor 
two  different  aspects  of  the  ionosphere's  vertical  electron 
concentration  profile,  Ne(h).  The  ionospheric  total  electron 
content  (TEC)  is  the  total  number  of  electrons  in  a vertical 
column  of  unit  cross  sectional  area  extending  through  the  iono- 
sphere. It  is  obtained  by  measuring  the  amount  of  Faraday 
rotation  a plane  polarized  wave  experiences  in  traveling  from  a 
satellite-borne  radio  beacon  to  a ground  receiving  station. 

The  peak  value  of  the  electron  concentration's  vertical  profile, 
denoted  Nmax>  ei.1  HmPg*;  is  obtained  from  ground-based  ionosonde 
measurements  of  the  ionosphere's  critical  (or  vertical  penetra- 
tion) frequencyt>ejfqr3T-  The  ratio  of  TEC/Nmax  is  called  the 
ionospheric  equivalent  slab  thickness , T".  a first-order  measure 
of  the  shape  of  the  Ne(h)  profile,  ^ 

An  introduction  to  the  relevant  parameters  and  the  physical 
principles  involved  (Chapter  1)  is  followed  by  a description  of 
a computer  model  for  ionospheric-plasmaspheric  structure  (Chap- 
ter II) . This  model  is  used  to  show  that  Faraday  rotation  can 
be  meaningfully  converted  to  TEC  at  L = 4 except  under  the  con- 
ditions of  a severe  trough,  a condition  which  is  primarily  a 
Winter  storm-time  phenomenon.  This  model  also  shows  that  gen- 
erally during  winter  nighttime  conditions  at  L = 4,  Faraday 
rotation  measures  the  total  electron  content  of  a vertical 
column  to  heights  greater  than  2000  km. 

Chapter  III  follows  with  a presentation  of  average  con- 
ditions at  L = 4,  obtained  from  over  three  years  of  Goose  Bay 
TEC,  St.  John’s  Nmax,  and  the  equivalent  slab  thickness  ob- 
tained via  T-  TEC/Nmax.  It  is  found  that  the  general  features 
of  the  ionosphere  are  the  same  as  those  found  at  mid  latitudes. 

In  Chapter  IV,  a particular  winter  nighttime  phenomenon 
near  L = 4,  increases  in  Nmax  accompanied  by  constant  or  de- 
clining TEC,  is  examined  in  detail.  A statistical  survey 
covering  24  years  of  December  foF2  data  and  two  October  to 
March  periods  reveals  that  nighttime  Nmax  increases  occur  very 
frequently  during  the  winter,  particularly  in  December.  Com- 
parable increases  in  TEC  appear  to  be  rare,  i.e.,  these  Nmax 
increases  are  most  often  accompanied  by  decreases  in  slab 
thickness.  It  is  shown  that  vertical  distortions  of  the 
ionospheric  Ne(h)  profile  are  the  major  cause  of  this  phenomenon, 
but  latitudinal  motions  of  the  electron  density  trough  may  also 
be  important  in  producing  this  effect.  These  vertical  and 
latitudinal  drift  motions  are  attributed  to  an  enhanced  magnet- 
ospheric  convection  electric  field. 


Unclassified 


SECU*lTy  CLASSIFICATION  OF  THIS  P AOEfNTien  Pmt»  F.nfrid) 


iii 


Acknowledgements 


We  are  grateful  to  the  Canadian  Department  of  Commu- 
nications for  providing  the  many  years  of  ionosonde  data 
from  St.  Johns,  Newfoundland,  and  the  ISIS-II  interpolated 
electron  density  topside  profiles,  used  in  this  study.  We 
also  wish  to  acknowledge  the  World  Data  Center-A  For  Solar 
Solar  Terrestrial  Physics  which  provided  multi-station  ground- 
based  ionosonde  values  and  the  true-height  analysis  necessary 
to  derive  the  bottom-side  electron  density  profiles  used  in 
our  case  studies.  Dr.  Jurgen  Buchau  of  the  Air  Force  Geo- 
physics Laboratory  most  graciously  performed  the  tedious 
scaling  of  several  ionograms  used  in  our  substorm  studies. 

Once  again,  John  Klobuchar  of  the  Air  Force  Geophysics 
Laboratory  contributed  insight  and  guidance  to  virtually 
every  aspect  of  the  investigations  summarized  in  this  report. 
His  continued  encouragment  in  these  matters  is  very  much 
appreciated. 


iv 


TABLE  OF  CONTENTS 

PAGE 

ACKNOWLEDGEMENTS  i;Li 

TABLE  OF  CONTENTS  iv 

LIST  OF  TABLES  vii 

LIST  OF  FIGURES  ix 

« 

CHAPTER  I THE  IONOSPHERE  AND  MAGNETOSPHERE  1 

1.1  The  Ionosphere  2 

1.1.1  Definition  of  the  Ionosphere  2 

1.1.2  Plasma  Frequency  and  Critical  Frequency  3 

1.1.3  How  the  Ionosphere  Is  Produced  and 

Maintained;  Ionospheric  Structure  4 

1.1.4  Ionospheric  Variations:  Trough,  Aurora  9 

1.2  The  Earth’s  Magnetic  Field  11 

1.3  The  Magnetosphere  13 

1.3.1  The  Magnetosphere:  Definition  and  Structure  13 

1.3.2  Relation  of  the  Ionospheric  Trough  to  the 

Plasmapause  17 

1.3.3  Geomagnetic  Storms,  Substorms  and  the 

Plasmapause  17 

CHAPTER  II  A MODEL  FOR  I ONOSPHER I C-PLASMASPHERI C 

STRUCTURE  19 

2.1  Description  of  a Model  for  lonospherlc- 

Plasmaspherl c Structure  20 


i 


2.1.1 

Faraday  Rotation 

20 

1.1.2 

The  Conversion  of  Faraday  Rotation  to  Total 

Electron  Content 

23 

2.1.3 

Model  Description 

25 

2. 1.3.1 

Model  Outline 

25 

2. 1.3. 2 

•'Exosphere1' 

27 

2. 1.3. 3 

"Ionosphere" 

29 

2.2 

Usefulness  of  Constant  M-factor  at  a High 

Latitude  Station  during  Winter  Night 

35 

2.  2. 1 

Ideal ized  Ar  rays 

35 

2.2.2 

Severe  Latitude  Gradient 

40 

2.2.3 

Conclusions 

47 

2.3 

lonosondes  and  Slab  Thickness 

48 

2.3.1 

lono  sondes 

48 

2.3.2 

Slab  Thickness 

50 

CHAPTER  III 

AVERAGE  BEHAVIOR  OF  THE  IONOSPHERE  NEAR 

SIXTY  DEGREES  MAGNETIC  LATITUDE 

52 

3.1 

Mean  and  Normalized  Monthly  Curves;  Median 

Behavi or 

53 

3.2 

Dl scusslon 

67 

CHAPTER  IV 

DISTORTIONS  OF  THE  WINTER  NIGHTTIME 

IONOSPHERE  NEAR  L=4 

73 

4.1 

Nighttime  Slab  Thickness  Decreases  Observed 

from  Goose  Bay  during  December#  1971 

74 

4.  2 

Statistical  Survey  of  Wi  nter  Nighttime 

Increases  In  foF2 

84 

Vi 


4.3 

Case  Studies:  December  9 and 

December  13,  1J71 

102 

4.3.1 

Description  of  Available  Data 

102 

4.  3.2 

Computer  Simulations  of  Slab  Thickness 

Decreases  for  Dec.  3 and  Dec.  15/  1971 

107 

4. 3.2. 1 

Preliminary  Considerations 

10  7 

4. 3. 2. 2 

Latitude  Shift  Method 

112 

4.3.2.  3 

Vertical  Drift  Method 

123 

4. 3. 2. 4 

Shift  and  Distort  Results 

141 

4.  3.  3 

Conclusions 

145 

4.4 

Mechanism  for  the  Slab  Thickness  Decreases 

14/ 

REFERENCES 


152 


vii 


J i i 

■ - 


LIST  OF  TABLES 


TABLE  PAGE 

11*1  Three  Faraday  Rotation  Observatories  35 

1 1-11  Uni  form  Latitude  Gradient  Input  Array  3b 

ll-lll  Best  h(R)  and  "%  dlf  at  420"  for  Uniform  Data 

and  Uniform  Gradient  Input  Arrays  37 

ll-IV  Severe  Latitude  Gradient  Input  Array  42 

II- V  Best  h(R)  and  ":£  dlf  at  420"  for  Severe 

Latitude  Gradient  Input  Array  45 

I I I-  I Data  Base  for  Mean  Monthly  Calculations  54 

I V-  I Data  Used  in  Statistical  Survey  05 

I V— 1 1 Statistical  Survey  of  Nocturnal  Winter  foF2 

Increases:  Seasonal  Analysis  53 


I V- I 1 1 Correlation  Coefficients  for  Seasonal  Analysis  95 
IV- IV  Statistical  Survey  of  Nocturnal  Winter  foF2 


Increases:  Solar  Cycle  Dependence  97 

IV-V  Correlation  Coefficients  for  Solar  Cycle 

Analysis  10U 

IV-VI  Calculated  and  Observed  TEC  110 

IV-VII  Observed  Increases-  December  9 113 

IV-VI  1 1 Observed  Increases-  December  15  113 

IV-IX  Latitude  Shift  Results-  December  9 

yzero*165kn  north  of  poleward  edee, 
yzero*215km  elesewhere  115 

IV-X  Latitude  Shift  Resu 1 ts-December  9 


I 


viii 

yzero*230km  north  of  poleward  edge* 
yzero»120km  in  trough  south  to  St.  John's* 
yzero*215km  from  St.  John's  southward 
IV-X1  Latitude  Shift  Results-  December  1j 

yzero»203km  north  of  poleward  edge* 
yzero=153km  In  trough  south  to  St.  John's* 
yzero*222km  from  St.  John's  southward 
IV-XII  Latitude  Shift  Resul ts-December  15 

yzeroa203km  north  of  poleward  edge* 
yzeroa12Jkm  In  trough  south  of  St.  John's, 
yzeroa222km  from  St.  John's  southward 
IV-XIII  Latitude  Shift  Results-  December  15 

yzero»230km  north  of  poleward  edge* 
yzero*120km  in  trough  south  to  St.  John’s* 
yzero*222km  from  St.  John's  southward 
IV-XIV  Summary  of  Changes  in  Nmax  and  fl  for 
Different  yzeros 

IV-XV  Vertical  Distortion  Results 

IV-XVI  Vertical  vs.  Equivalent  Vertical  TEC- 

Distortion  Method 

IV-XVII  Shift  and  Distort  Results 

IV-XVI I I Magnetospher i c Dawn-Dusk  Electric  Field 


Values  for  Varying  Levels  of  Geomagnetic 
Activity*  from  (Carpenter*  1970) 


150 


ix 


LIST  OF  FIGURES 

FIGURE  PAGE 

1-1  A typical  daytime  ionospheric  profile 

from  ( Papagi annl s,  1172)  8 

1-2  Magnetic  dipole  geometry  12 

1-3  Principal  features  of  the  magnetosphere 

from  (Heikkl la,  1973)  14 

1-4  Equatorial  radius  of  the  plasmapause  vs.  local 
time  during  periods  of  moderate,  steady 
geomagnetic  agitation  1j 

I- 5  Nishida's  (196G)  magnetospher  i c convection 

pattern  16 

I 1-1  Relation  between  distance  along  ray  path 

and  height  23 

II- 2  A sample  page  of  Isis  II  Interpolated  Ne(h)  data  30 

11-3  Geometry  showing  relation  between  station  latitude 

and  height  at  which  the  ray  path  crosses  the 
lowest  invariant  latitude  39 

11-4  Distances  along  ray  path  where  8U2  of  Q is 

accumulated  (10 'i  above,  10%  below)  for  case  of 
uniform  gradient  41 

11-5  Isis  II  pass,  Dec.  18,  1971  43 

II- 6  Electron  density  vs.  height  along  ray  path; 

Dec.  18,  1971;  +8*  degree  shift  46 

II I- 1  Mean  monthly  TEC  diurnal  curves 


55 


X 


I I I -2  Mean  monthly  Nmax  diurnal  curves  56 

MI-3  Mean  monthly  Slab  Thickness  diurnal  curves  57 

II I- 4  Normalized  mean  monthly  TEC  diurnal  curves  uO 

I I 1-5  Normalized  mean  monthly  Nmax  diurnal  curves  61 

MI-6  Normalized  mean  monthly  slab  thickness 

diurnal  curves  62 

MI-7  Maximum  and  minimum  mean  monthly  TEC  u4 

MI-8  Maximum  and  Minimum  mean  monthly  normalized  TEC  65 

II  I -9  Median  monthly  TEC  contour  plot  u6 

IV- 1  Nmax  and  TEC  for  five  nights  when  the  increasing 

Nmax-cons tant  TEC  effect  occurred  75 

IV-2  Several  parameters  plotted  for  the  night  of 

Dec.  8,  1971  77 

I V- 3 Several  parameters  plotted  for  the  night  of 

Dec.  9,  1971  78 

IV-4  Several  parameters  plotted  for  the  night  of 

Dec.  10,  19  7i  79 

IV-5  Several  parameters  plotted  for  the  night  of 


Dec.  12,  19/1  80 
I V— 6 Several  parameters  plotted  for  the  night  of 

Dec.  15,  19/1  81 
I V—  7 Overplots  of  Nmax  and  TEC  for  December,  1971  82 
IV-8  Classification  example:  Low  bump,  plateau  68 
I V—  9 Classification  example:  bump,  short  bump  89 
IV-10  Classification  example:  large  bump, i rregu lar  bump  30 
I V—  1 1 Classification  example:  no  bump  91 


wmmmrnmm 


X I 

I V— 12  Classification  example:  Insufficient  data  92 

IV-13  Statistical  survey,  seasonal  analysis,  fractional 

occurrence  of  all  foF2  increases  94 

IV-14  Statistical  survey,  solar  cycle  dependence, 

fractional  occurrence  of  all  foF2  increases  99 

IV-15  lonograms  taken  during  the  Dec.  15  event  104 

IV-lb  Ne(h)  vs.  latitude:  Isis  II  pass  4:20UT,  Dec.  9 106 

IV-17  Ne(h)  vs.  latitude:  Isis  II  pass  2:25UT,  Dec.  15  108 

IV-18  Model  Nmax  and  Cl  for  Dec.  15  using  various 

yzeros  for  a wide  range  of  latitude  shifts.  120 

IV-19  We(h)  vs.  latitude  contours  for  Dec.  13,  0 and 
-6  degree  shifts,  using  yzeros  shown  in 
Table  I V-XI  I I 121 

IV-20  Ne(h)  vs.  latitude:  Isis  II  pass  b:laUT,  Dec.  9 124 

IV-21  Ne(h)  vs.  latitude:  Isis  II  pass  4:20111,  Dec.  15  125 

IV-22  Typical  nighttime  Ne(h)  profile  showing  that 

N(hmax-Aji)  < N(hmax+Zih)  150 

IV-23  Ne(h)  vs.  latitude:  Dec.  9 distorted  array  155 

IV-2h  Ne(h)  vs.  latitude:  Dec.  15  distorted  array  13b 

I V— 2 5 Ne(h)  along  ray  path:  Dec.  9 original  and 

distorted  array  13a 

IV-26  Ne(h)  along  ray  path:  Dec.  15  original  and 

distorted  array  139 

IV-27  Hypothetical  response  of  Nmax  to  a vertical 
distortion  at  two  different  times  during 
the  night  14  7 


2 


1.1  The  Ionosphere 
1,1.1  Definition  of  the  Ionosphere 

The  Ionosphere  Is  that  part  of  the  Earth's  atmosphere 
which  contains  a sufficient  density  of  Ionized  particles  to 
affect  the  propagation  of  radio  waves.  Based  on  this 
definition  the  Ionosphere  can  be  said  to  extend  from  a 
height  of  about  40  km  to  several  earth  radii  (Davies*  1965) 
although  the  term  Ionosphere  is  not  commonly  used  in 
reference  to  heights  above  a few  thousand  kilometers.  This 
definition  emphasizes  the  practical  nature  of  ionospheric 
research.  The  density  of  Ionization  Is  only  about  one 
thousandth  that  of  the  neutral  particles  yet  this  relatively 
Insignificant  component  of  the  Earth's  upper  atmosphere  has 
made  possible  long  distance  radio  communication  and  the 
variations  In  Its  behavior  are  thus  of  Interest  to  a wide 
range  of  people.  The 
Ionosphere  not  only  has  great  practical  significance  but  is 
of  Interest  also  because  of  Its  scientific  research  value  as 
a conveniently  located  cosmic  plasma  laboratory  for  diverse 
physical  processes.  Because  of  the  effeet  on  radio  waves  "f 
the  Ionospheric  medium  through  which  they  propagate*  the 
most  common  methods  for  observing  the  Ionosphere  use  radio 
waves*  and  the  most  frequently  used  radio  technique  is  that 
of  vertical  sounding  (see  sec,  2,3.1). 
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1.1.2  Plasma  Frequency  and  Critical  Frequency 

The  effect  of  Ionospheric  plasma  on  radio  waves  Is  due 
to  the  electromagnetic  nature  of  the  radiation.  The  electric 
field  of  the  passing  radio  wave  sets  the  light  electrons 
oscillating  at  the  frequency  of  the  wave.  The  electrons 
absorb  energy  from  the  wave  and  reradlate  It.  The  electrons 
thus  act  as  forced  oscillators,  and  like  a block  on  a spring 
they  have  a natural  frequency  of  vibration  known  as  the 
plasma  frequency,  fN , It  Is  easy  to  show  that: 

f» * 3 ^ ai) 

where  Is  In  megahertz  and  N Is  the  electron  density 
measured  In  units  of  10*e1/cm*  ( Papaglannt s,  1972).  Thus  the 
plasma  frequency  depends  only  on  the  electron  density.  If  a 
radio  wave  Is  Incident  vertically  at  or  below  this  plasma 
frequency  the  wave  Is  totally  reflected  while  at  frequencies 
much  larger  than  the  plasma  frequency  the  Incident  wave  Is 
transmitted  with  little  attenuation. 

The  concentration  of  electrons  throughout  the 
Ionosphere  Is  far  from  uniform.  It  Is  negligible  at  the 
surface  of  the  earth,  also  very  small  a few  earth  radii 
away,  and  reaches  a maximum  normally  between  200km  and  500km 
In  height.  This  height  Is  given  the  symbol  hmax  (or  hmF2), 
while  the  electron  concentration  at  hmax  Is  called  the  peak 
density  and  given  the  symbol  Nmax  or  NmF2  (see  section 
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1,1,3).  This  peak  density  corresponds  to  a “critical*'  plasma 
frequency  denoted  foF2  (from  equation  1,1),  which  Is  the 
maximum  frequency  for  a vertically  Inclder  radio  wave  to  be 
reflected.  The  measurement  of  critical  frequency  through  the 
use  of  lonosondes  (see  section  2.3,1)  directly  provides 
Nmax,  Because  of  this  foF2  Is  perhaps  the  most  widely 
studied  Ionospheric  parameter. 

1.1.3  How  the  Ionosphere  Is  Produced  and 
Maintained;  Ionospheric  Structure 

The  main  constituents  of  the  neutral  atmosphere  'are 
atomic  and  molecular  oxygen  and  nitrogen,  atomic  helium  and 
hydrogen.  Helium  Is  not  an  Important  constituent  below  about 
500km  while  hydrogen  Is  unimportant  below  about  1000km. 
Oxygen  begins  to  be  dissociated  above  100km  by  ultraviolet 
radiation.  It  then  diffuses  upward  making  It  the  dominant 
source  of  Ions  In  the  F-reglon  (see  below).  Molecular 
nitrogen  Is  not  readily  dissociated  but  atomic  nitrogen  as 
well  as  nitric  oxide  Is  produced  by  the  Ion-atom  Interchange 
and  dissociative  recombination  reactions  described  below. 

In  general,  tons  and  electrons  are  created  from  neutral 
particles  through  the  action  of  ultraviolet  and  x-ray 
radiation,  with  a small  contribution  from  cosmic  rays.  The 
rate  of  production  of  Ions  and  electrons  In  the  Earth's 

upper  atmosphere  Is  proportional  to  the  number  of  neutral 

I 


particles  available  to  be  Ionized  and  also  to  the  Intensity 
of  the  Ionizing  radiation.  At  the  greatest  heights  there  Is 
very  little  gas  to  be  Ionized  while  below  a height  of  about 
100km  most  of  the  I net  dent  radiation  has  already  been 


absorbed.  Thus  the  maximum  production  rate  will  occur  where 


neither  the  radiation  nor  the  concentration  of  neutral 


particles  Is  vanishingly  small 


Loss  of  Ionization  by  recombination  with  Ions  through  a 


reaction  of  the  form 


where  A*  Is  an  atomic  ton,  will  not  conserve  both  energy  and 
momentum  except  In  the  presence  of  a neutral  third  particle. 


The  concentrations  of  neutral  particles  are  sufficiently 


large  In  the  lower  Ionosphere  (D-reglon,  E-reglon,  and 


Fl-reglons)  so  that  this  direct  recombination  process  Is  the 


dominant  one 


The  rate  of  recombination  depends  on  Ion  and  electron 


concentrations,  which  are  equal,  and  so  it  depends  on  the 


square  of  the  electron  concentration.  Thus 


where  at  Is  the  recombination  coefficient 


Above  the  Fl-reglon  lies  the  F2-reglon  where  the 


density  of  neutral  particles  has  fallen  to  very  low  levels 


so  that  three  body  collisions  become  extremely  rare.  Then 


more  complicated  loss  processes  Involving  combinations  of 


two  body  collisions  between  neutral  particles.  Ions 
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electrons,  become  most  Important.  Notice  that  the  reaction: 

e"  -r  AB*  A + B U-'V) 

will  conserve  both  energy  and  momentum.  This  process  Is 
known  as  dissociative  recombination.  Atomic  oxygen  Ions  are 
often  recombined  with  electrons  by  first  undergoing  Ion-atom 
Interchange  followed  by  dissociative  recombination  as 
fol I ows : 

0*  + Ni  -»  N0+  + N (1.5) 
e-  * NO*—*  N f 0 (l.i) 


It  can  also  first  undergo  charge  transfer  with 
molecular  oxygen  followed  again  by  dissociative 
recomb  1 nation: 


0+  ♦ o^— >•  0 + 0+  (1.7) 

e-  + 0*-*  0*0  (1.8) 


Equations  (1,5)  and  (1,7)  proceed  more  slowly  than 
(1.6)  and  (1.8),  respect  I vel y,  because  the  concentrations  of 
neutral  particles  are  so  low.  Therefore,  the  upper 
Ionospheric  loss  processes  are  dependent  upon  the 
concentrations  of  the  neutrals.  For  these  loss  processes 


then : 


dU  . _ 

dt 


(: n ) 


where  p Is  called  the  attachment  coefficient.  Equation  (1.9) 
Is  In  the  form  of  a simple  “at tachment- 1 1 ke"  process  as  If 
electrons  were  lost  by  simply  attaching  themselves  to  the 
rare  neutral  particles. 
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The  continuity  equation  describing  N ( t ) Is: 

= q _L  0.10) 

ot 

where  q Is  the  production  rate  and  l the  loss  rate.  In  his 
classic  paper  (Chapman,  1931),  Sidney  Chapman  examined  the 
steady  state  case  dN/dt*0,  whence  q»L.  By  making  several 
other  assumptions  he  showed  that  for 

S * L = - 

N * M**  expill  - * - secX«p(-ifl  0.1 1) 

while  for  q * L » 

N-  2 ' ietXwp('^] 

where  X * zenith  angle,  z - (h-hmax)/H  with  H ■ a constant 
scale  height. 


U.ii) 


The  assumptions  made,  viz,,  an  atmosphere  with  a single 
constituent,  plane  stratification,  a parallel  beam  of 
monochromatic  solar  Ionizing  radiation  and  an  Isothermal 
atmosphere,  limit  th<*  usefulness  of  Chapman's  theory,  but 
Its  value  Is  evident  from  Its  continued  widespread  use 
particularly  In  simple  Ionospheric  models. 

The  various  regions  of  the  Ionosphere  seen  during  the 
daytime  are  depicted  In  figure  1-1. 


(■ 


[ 


•: 


figure  1-1  A typical  daytime  Ionospheric  profile 
from  (Papaglannts#  1972). 


While  the  production  rate  reaches  Its  maximum  In  the 
Fl-reglon#  the  maximum  concentration  of  electrons  is  found 
In  the  F2-reglon  where  the  loss  rate  Is  so  small  that  some 
plasma  remains  even  at  night.  Distinct  D#  E#  and  FI  regions# 
however#  disappear  at  night.  While  the  loss  rate  in  the 
F2-reglon  Is  low  one  would  expect  from  (1,9)  a steady  decay 
In  the  course  of  a night.  While  the  minimum  value  of  Nmax  is 
usually  found  just  before  dawn#  there  Is  often  not  a smooth 
nighttime  decay.  The  problem  of  explaining  the  maintenance 
of  the  nighttime  Ionosphere  has  been  Investigated  for  many 
years.  One  aspect  of  the  problem-  certain  winter  nighttime 
Increases  In  foF2-  will  be  taken  up  In  Chapter  4, 
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1,1.4  Ionospheric  variations:  Trough,  Aurora 

Many  Ionospheric  variations  and  "anomalies"  (so  called 
because  they  are  not  predicted  by  simple  Chapman  theory) 
have  been  studied.  It  would  not  be  appropriate  to  even 
briefly  describe  all  these  variations  (see,  for  Instance 
(Papagtannls,  1972),  section  2.6),  but  a few  of  greatest 
Interest  to  this  study  should  be  mentioned. 

The  seasonal  anomaly  and  December  anomaly  are  the 
higher  values  of  noontime  foF2  In  winter  than  in  summer  and 
the  generally  higher  foF2  In  December  than  In  other  months, 
respectively. 

The  latitudinal  dependence  of  the  electron  density  Is 
primarily  due  to  varying  elevation  of  the  noontime  sun  with 
latitude.  The  midlatitude  trough  has  a more  complicated 
Interpretation,  however,  which  will  be  discussed  In  section 
1.3.2, 

The  11  year  solar  cycle  Is  a regular  variation  In  solar 
activity  as  well  as  In  the  average  intensity  of  UV  and 
x-radlatlon  which  Is  accompanied  by  a corresponding 
variation  in  Ionospheric  plasma  density. 

The  disturbed  Ionosphere  Is  character  I zed  by  such 
phenomena  as : 

Sporadlc-E,  which  Is  the  formation  of  a thin  layer  of 
Ionization  between  100  and  140  km  In  height  (Pike,  Wagner 
and  Akasofu,  1975), 
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Spread-F,  the  large  spread  In  height  from  which 
lonosonde  echoes  are  returned  due  to  a MblobbyM  F-reglon. 

Ionospheric  storms,  significant  enhancements  or 
depletions  In  Ionospheric  total  electron  content  and  foF2 
often  associated  with  geomagnetic  storms  (see  section 
1,3.3). 

Aurora  cause  significant  Increases  In  electron  density 
parti culary  In  the  E-reglon.  While  these  phenomena  occur  at 
high  latitudes,  they  are  often  found  at  ml d-1 atl tudes  during 
occasional  Intense  magnetic  storms. 


1,2  The  Earth's  Magnetic  Field 


For  simplicity,  the  Earth's  magnetic  field  Is  usually 
approximated  by  a centered  dipole.  This  approximation  Is 
nearly  perfect  at  distances  of  a few  Earth  radii,  usually 
gives  a value  within  1 % of  the  real  field  near  the  Earth's 
surface  (except  during  a great  magnetic  storm  when  the 
dlscrepency  may  be  3£),  but  becomes  poor  beyond  about  8 
earth  radii  or  L=8  (the  equatorial  distance  expressed  in 
Earth  radii  Is  called  the  L-value), 

Some  of  the  basic  equations  describing  this  centered 
dipole  magnetic  field  follow, 

B - -$-0  + 3 sin'e)*  (1.13) 

where  B Is  expressed  In  gauss,  0 is  the  magnetic  latitude,  r 
the  equatorial  distance,  and  M“8,05(  + .002)xl0*'5gauss,  Also: 

tan  I ® £tan0  (1.14) 

where  I Is  the  Inclination  of  the  field  line  with  respect  to 
the  vertical.  One  characteristic  of  magnetic  field  lines  is 
that  the  Inclination  remains  constant  for  all  field  lines  at 
a given  magnetic  latitude. 

For  a given  field  line: 

r-  r.  cot*0  s LReCos'9  (I.J5) 

where  r0  Is  the  distance  to  the  intersection  of  the  field 
line  with  the  magnetic  equator,  and  Re  is  the  Earth's 


radius.  For  r«Re: 


where^ls  called  the  Invariant  latitude. 

Since  r ■ Re  ♦ h : 

cosxA'  = L+  (j.17; 

wherey^  Is  called  the  generalized  Invariant  latitude.  For 
example*  at  h*1000km  and  L=4*  «57°28'.  Figure  1-2  shows 

some  of  the  magnetic  dipole  geometry. 


figure  1-2  Magnetic  dipole  geometry 
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1,3  The  Magnetosphere 


1,3,1  The  Magnetosphere:  Definition  and  Structure 


The  magnetosphere  has  been  defined  as:  "...the  region 
above  the  Ionosphere  In  which  the  magnetic  field  of  the 
earth  has  a dominant  control  over  the  motions  of  gas  and 
fast  charged  particles*'  (Gold,  1959),  The  lower  boundary  of 
the  magnetosphere,  l,e,,  Its  boundary  with  the  Ionosphere, 
Is  traditionally  near  1000km  (Parker,  1968).  Even  In  the 
F-reglon  (near  300km),  however,  the  presence  of  the  magnetic 
field  exerts  a strong  Influence  on  the  plasma,  re:  the 
"frozen  field"  concept  (Rlshbeth  and  Garrlott,  1969,  page 
244-5),  and  as  we  shall  see  In  Chapter  IV,  the  magnetic 
field  can  cause  substantial  drift  motions  of  the  Ionospheric 
plasma.  In  1959  the  magnetosphere  was  thought  to  extend 
about  10  Earth  radii  In  the  equatorial  plane  symmetrically 
with  respect  to  local  time.  A more  recent  diagram  showing 
the  remarkable  progress  In  the  scientific  Investigation  of 
the  magnetosphere  which  has  paralleled  the  advent  of  the 
space  age  Is  given  In  figure  1-3, 

It  was  the  effect  of  the  solar  wind  on  the  near  earth 
environment  which  most  radically  altered  the  original 
picture  of  the  magnetosphere.  The  pressure  of  this  solar 
wind  leads  to  a compression  of  the  dayslde  magnetic  field 
lines  and  the  long  streaming  open  field  lines  of  the 


geomagnetic  tall. 

The  boundary  on  the  nlghtside  between  field  lines  which 
are  always  closed  and  those  which  are  sometimes  open  has 
been  associated  with  the  equatorial  piasmapause,  discovered 
tn  1963  (Carpenter,  1963).  Usl ng  whistler  data.  Carpenter 
found  a sudden  decrease  from  about  lOOOel/cm*  to  only  about 
lei /cm*  at  this  plasma  knee  or  piasmapause.  Carpenter 
studied  the  temporal  dependence  of  piasmapause  position  and 
found  that  near  midnight  It  Is  located  at  L“4  (Carpenter, 
1966).  In  the  course  of  a day  the  piasmapause  position  was 
seen  to  vary  from  a little  over  L=3  to  over  L»5  (see  figure 
f-4). 

IX 


figure  1-4  Equatorial  radius  of  the  piasmapause 
vs.  local  time  during  periods  of  moderate  steady 
geomagnetic  agitation;  (Carpenter,  1966). 

The  variations  In  the  position  of  the  piasmapause  have 
been  dealt  with  In  the  context  of  theories  of  magnetospher I c 
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structure.  Even  before  Carpenter's  discovery  of  the  plasma 
"knee",  Axford  and  Hines  (1961)  suggested  that  sunward 
magnetospherlc  convection  would  be  driven  by  a dawn  to  dusk 
electric  field  (E»»)  caused  by  the  interaction  of  the  sotar 
wind  with  the  magnetosphere,  Nlshlda  (1966)  showed  that  a 
convection  pattern  would  develop  with  outward  moving  plasma 
obtained  from  the  draining  of  magnetic  tubes  of  force  which 
open  as  they  drift  across  the  magnetic  tali,  and  Inward 
moving  plasma  obtained  from  the  ExB  mechanism.  The  boundary 
between  tubes  of  force  which  open  and  those  which  never  do 
Nlshlda  Identified  with  the  plasmapause.  His  convection 
pattern  Is  reproduced  In  figure  1-5,  Brice  (1967)  further 
developed  the  theory  of  macnetospher l c convection  by 
examining  It  from  the  point  of  view  of  field  line  merging. 


SOLAR  WIND 


figure  1-5 


from  (Nlshlda, 1966) 
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1.3.2  Relation  of  the  Ionospheric  Trough  to 
the  Plasmapause 

Since  charged  particles  are  generally  “tied  to" 

magnetic  field  lines,  the  field  line  Intersecting  the 
equatorial  plasmapause  may  be  considered  also  to  be  the 
three-dimensional  boundary  of  the  plasmasphere,  a region  of 
much  higher  electron  density  than  the  outer  magnetosphere. 
This  same  field  line  Intersects  the  surface  of  the  Earth  at 
the  Invariant  latitude  given  by  (1,1(3).  One  then  expects  the 
Ionospheric  plasma  density  to  drop  sharply  at  that  latitude. 
In  fact  deep  Ionospheric  troughs  are  found  at  or  near 
plasmapause  L-values,  Electron  density  may  drop  by  a factor 
of  two  within  10km  at  a height  of  1000km,  (see  (Rlshbeth  and 
Garriott,  1969),  sections  5,52  and  7.54).  The  term 

plasmapause  no  longer  refers  merely  to  the  equatorial 
"knee,"  but  Is  now  generally  used  In  reference  to  the  field 
line  which  serves  as  the  plasmasphere  boundary, 

1.3.3  Geomagnetic  Storms,  Substorms  and 
the  Plasmapause 

The  general  features  of  a geomagnetic  storm  may  be 
outlined  as  follows, 

A solar  flare  ejects  a large  cloud  of  plasma  which  In  a 
few  days  Impinges  on  the  Earth's  magnetosphere.  The 


resultant  sudden  compression  of  the  geomagnetic  field  lines 
usually  causes  a sudden  Increase  In  the  observed  magnetic 
field  strength  known  as  a sudden  storm  commencement.  The 
Increase  In  magnetic  field  strength  Is  referred  to  as  the 
positive  phase  of  the  storm.  Enhanced  magnetospher I c 
convection  leads  to  a buildup  of  particles  at  about  L*4 
known  as  a ring  current  which  depresses  the  magnetic  field 
by  hundreds  of  gammas  (1  gauss*l00000  gammas).  This 
depression  is  referred  to  as  the  storm's  main  phase.  There 
Is  then  a gradual  recovery  which  lasts  2 to  3 days. 

A substorm  has  no  positive  phase  and  is  distinguished 
by  a short  main  phase  (1-3  hours).  It  is  believed  that  the 
main  phase  of  a geomagnetic  storm  is  actually  the 
superposition  of  several  substorms  (Akasofu,  1968). 

During  magnetically  disturbed  periods,  the  plasmapause 
is  seen  to  move  earthward  accompanied  by  a southward  motion 
of  the  ionospheric  trough.  Carpenter  (1966)  found  about  a 
six  hour  phase  lag  between  Kp  (a  popular  index  of 
geomagnetic  activity)  and  inward  plasmapause  motion.  Brace 
et.  al.  (1974)  noted  that  during  the  August  1972  geomagnetic 
storm  the  plasmapause  moved  earthward  from  1-4.5  to  L-1.9  at 
the  peak  of  the  storm,  Tulunay  and  Sayers  (1971)  completed  a 
statistical  analysis  of  over  1000  troughs  which  showed  an 
inverse  correlation  between  Kp  and  trough  latitude* 
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2,1  Description  of  a Model 
for  lonospher I c-PI asmaspher l c Structure 


2.1,1  Faraday  Rotation 


The  Index  of  refraction  in  the  presence  of  a magnetic 
field  and  ionized  plasma  with  low  collision  frequency  is 
found  from  the  fundamental  Appleton-Har tree  equation  and  Is 
given  by: 


i _ xii-x) 

1 (l-X)-i'C 


u.i) 


where  X * fN*  /f4  (2.2) 

and  Yt  * YslnB,  YL  =Ycos8  (2,3) 
where  Y * fH  /f  (2.4) 

Here  f * wave  frequency  and 

f . eH 

w W air  me 

* R 

•gyrof requency  ~ — (2.5) 

Jtirmc 

and  8 ■ propagation  angle  • angle  between  the  magnetic  field 
and  the  wave  normal.  Unless  8 is  close  to  9d*  a considerable 
simplification  can  be  made.  We  can  take 

*rT4  « 


whence: 


This  Is  known  as 


= I 


- _x 


U.fc) 


I tXu 

the  "quasi - longl tudi nal 


approximation."  1 1 


is  excellent  for  8<86  at  f»100MHz  for  Instance 
(Haj  eb-Hossei  ni  eh#  19  73).  At  lower  frequencies  It  Is  a 
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poorer  approximation#  of  course#  but  Faraday  rotation 
studies  are  normally  carried  out  using  frequencies  larger 
than  100MHz  so  Its  use  Is  generally  justified. 

The  magnetosphere  is  a so-called  “magnetolonl c“  medium. 
Such  a medium#  like  a crystal#  has  blrefringent  properties# 
l.e.#  a radio  wave  Is  broken  up  Into  two  components#  an 
ordinary  wave#  and  an  extraordinary  wave.  The  Faraday 
rotation#  a rotation  of  the  plane  of  polarization  of  a radio 
wave#  is  best  explained  by  using  the  fact  that  a linearly 
polarized  wave  can  be  regarded  as  the  sum  of  two  circularly 
polarized  waves  rotating  in  opposite  directions#  one  of 
which  may  be  Identified  with  the  ordinary  wave#  the  other 
with  the  extraordi nary  wave.  The  ordinary  and  extraordl nary 
waves  have  different  indices  of  refraction  corresponding  to 
the  positive  and  negative  signs  respectively  in  C2.1)  and 
(2.6).  We  also  have: 

5 C/VpK  U.1) 

where  Is  the  phase  velocity.  Thus  the  phase  velocities  of 
the  two  modes  differ  resulting  In  a phase  lag: 

an  * iu*.  - <h.)  a.© 

where: 

^ • af  ■ 'w 

J[i  - SUL  (r.  (wo) 

c 


Thus : 
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But  (2*6)  can  be  approximated  as  follows: 

/»*-  V - iTyf  * I -X(UYt) 

/»  = [(-  XO  * Y0lfc  « 1-4X11*  Yu)  an) 


Then: 

/"•«»-  iXll-Yu) 

/*«*!-  iXll+Yu) 


Substitution  Into  (2.10)  gives: 

an=ifiixY.  u.i» 

C 1 

Using  (2. 2), (2. 3), (2. 5),  and  (1.1)  we  arrive  at: 


<m  * 


A (in) 

here  K • .236  In  CGS  units  with  B expressed  in  gammas. 

From  (2.13)  one  can  see  that  knowledge  of  the  Faraday 
rotation  H ^ * radio  wave  transmitted  from  a satellite 
provides  information  about  the  total  amount  of  plasma 
encountered  along  the  ray  path.  It  would  be  most  useful  to 
translate  the  Information  contained  in  (2.13)  into  knowledge 
of  a vertical  column  of  standard  cross  section  located  at  a 
certain  point  along  a slanted  ray  path.  The  number  of 
electrons  In  such  a column  Is  referred  to  as  the  total 
electron  content  of  the  column  (TEC).  We  will  next  outline 
the  technique  for  deriving  the  TEC  from  O. 
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2.1.2  The  Conversion  of  Faraday  Rotation  to 
Total  Electron  Content 


n 


= -^(VweWcosesttXA 


f*). 


- * (kMNWJh 
f 


(1.14) 


where: 


M-  B(W  cos  8 setX  (**5) 

For  a geostationary  satellite  (at  r»6.619Re)#  the  upper 
limit  of  the  Integral  Is  usually  taken  to  be  1000km  or 
2000km  because  of  the  small  amount  of  Faraday  rotation  which 
occurs  In  the  plasmasphere  and  beyond  (Tltherldge,  1972). 
Over  the  small  range  from  the  surface  of  the  earth  to  1000 
or  2000km,  M changes  slowly  so  a mean  value  Is  taken 
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corresponding  to  the  height  along  the  ray  path  at  which  the 
vertical  column  is  taken.  Then: 

SI  * -k  R (kN(U<lh  ( 3i.lt) 

f Jo 

The  integral  In  (2,16)  is  referred  to  as  the  equivalent 
vertical  TEC,  In  the  absence  of  severe  latitude  gradients 
the  equivalent  vertical  TEC  would  be  expected  to  be  quite 
close  to  the  actual  total  electron  content  of  the  vertical 
column. 

A choice  must  be  made  of  the  height  along  the  ray  path 
h CM)  at  which  the  vertical  column  will  be  placed  This 
choice  is  made  on  the  basis  of  model  studies  (such  as  that 
described  in  the  next  section).  The  single  height  is  chosen 
at  which  the  equivalent  vertical  TEC  and  the  actual  electron 
content  of  the  vertical  column  are  in  best  agreement  for  a 
wide  variety  of  conditions,  i.e,,  solar  cycle,  seasonal, 
diurnal  and  latitudinal  variations,  for  instance,  must  be 
examined  by  such  a model. 

Using  such  a model  Titheridge  (1972)  concluded  that  use 
of  h(R)»420km  with  Faraday  rotation  measurements  of 
geostationary  satellite  signals  (satellite  at  r"6,619Re  in 
the  equatorial  plane)  would  give  the  electron  content  up  to 
2000km  of  a vertical  column  to  5%  accuracy  for  a wide 
variety  of  conditions.  Little  Faraday  rotation  is  added 
above  2000km.  If  one  wished  to  extend  the  vertical  column  to 
the  height  of  the  satellite,  no  single  value  of  h(R)  would 
suffice.  This  results  from  the  high  variability  In  electron 
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density  In  these  upper  regions  due  to  the  wide  variations  in 
the  height  of  the  0+/H+  transition  level.  This  can  lead  to 
changes  in  by  even  a factor  of  four,  corresponding  to  a , 
3000km  variation  In  h(M).  Papaglannis  et.  al.  (1975)  showed 
that  the  M value  used  to  calculate  the  electron  content  up 
to  1000km  from  signals  transmitted  by  a satellite  at  1000km 
is  approximately  the  same  (within  l-2£)  as  the  M value  used 
to  calculate  the  electron  content  up  to  2000km  from  signals 
transmitted  by  a geostationary  satellite. 

2.1.3  Model  Description 
2. 1.3.1  Model  Outline 

The  computer  model  used  in  this  study  divides  the 
Earth's  environment  into  two  regions:  the  regions  below  and 

above  1000km.  In  the  lower  region  which  may  be  loosely 
called  the  "Ionosphere"  we  create  a two  dimensional  "real 
data"  array  by  reading  into  the  program  data  obtained  from 
the  Isis  II  topside  sounder  satellite,  i nterpolat I ng  between 
the  data  points  and  adding  a model  bottoms! de  profile.  For 
the  upper  region  (the  "exosphere")  we  can  calculate  the 
electron  density  at  any  location  by  using  one  of  two  models, 
each  dependent  upon  conditions  at  a reference  height 
(1000km).  The  choice  between  the  two  models  : a 

"collisionless"  (r-*  ) model  or  a diffusive  equilibrium 


model#  depends  upon  whether  the  point  in  question  lies 
within  or  is  outside  the  pi asmasphere.  The  Faraday  rotation 
and  TEC  accumulated  along  a ray  path  down  to  lOUOkm  In 
height  is  added  up#  and  combined  with  that  of  the  lower 
region.  The  Faraday  rotation  for  the  lower  region  is 
obtained  by  interpolating  (or  extrapolating)  between  the 
“real  data  arr ay'1  points.  So  by  matching  the  two  sections 
one  can  follow  the  ray  path  from  the  satellite  to  the 
ground.  The  vertical  TEC  below  1000km  at  various  points 
along  the  ray  path  Is  calculated  by  the  creation  of  a 51  by 
51  Mi  nterpolated1'  density  array  which  divides  the  range  of 
latitude  the  ray  path  covers  below  1000 km  into  51  vertical 
segments#  one  vertical  column  through  each  of  51  height 
steps  of  20km. 

A more  detailed  description  follows. 
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2. 1.3. 2 “Exosphere" 

Anger am I and  Thomas  (1964)  studied  theoretically  the 
distribution  of  pi asmaspher I c electrons  and  ions.  For  the 
case  of  a constant  temperature  along  a given  field  line  and 
equal  ion  and  electron  temperatures,  their  results  give: 


N(s)  = [ exp(-i/Hl)  4 Vxp(-2/4Hl)  »1,exp(-: 


N 


eo 


i 


1 


an) 


Here  Neo  * the  electron  density  at  a reference  height 
(1000km  in  this  study)  along  the  field  line.  N(s)  is  the 
electron  density  at  a given  distance  along  the  field  line. 
HI  is  a constant  scale  height.  where  the  etas  are 
the  ratios  of  the  ion  densities  hydrogen,  helium  and  oxygen 
to  that  of  oxygen.  The  parameter  z,  called  the  geopotential 
height.  Is  given  by: 

* - MO -'£&]♦  u-,l) 

where  R0  Is  the  distance  from  the  center  of  the  Earth  to  the 
reference  height,  0O  Is  the  generalized  invariant  latitude 
at  the  reference  height,  0'  is  the  magnetic  latitude  of  the 
given  point  along  the  field  line,  W Is  the  Earth's  angular 
velocity  of  rotation  and  g#  Is  the  gravitational 
acceleration  at  the  reference  level. 

In  utilizing  these  equations,  the  temperature  (needed 
for  the  scale  height)  and  the  reference  level  concentrations 
of  electrons  and  relative  Ion  concentrations  must  be 
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provided.  The  electron  density  at  our  chosen  reference  level 
(1000km)  is  obtained  by  interpolating  linearly  (or 
extrapolating)  between  the  Isis  ii  real  data  points.  For  the 
other  four  parameters  the  following  distributions  were 
assumed  for  the  winter  nighttime  conditions  we  will  be 
considering  later  (the  ion  concentrations  are  fractions): 
temp  « 1125.0  ♦ 1050. 0(A'  - 30#.Q)/3U°.0  (2.19) 

oxygen  - 0.10  ♦ 0.80(A'  * 3**.0)/22#.0 
hydrogen  - 0.  82  - 0.  73(A'  - 33,.0)/2/.0  (2.20) 

helium  • 1 - oxygen  - hydrogen 
except  that  for  A* (which  refers  to  1000km)  greater  than  55*: 
oxygen  • 0.90 

hydrogen  ■ .09  (2.21) 

helium  * .01 

The  temperature  relation  was  scaled  from  (Brace,  Reddy, 
and  Mayr,  1967),  while  the  composition  formulas  (2.20)  and 
(2.21)  were  obtained,  following  (Mendillo  and  Klobuchar, 

1974) ,  from  Millstone  Kill  data  and  from  the  observations  of 
Prasad  (1970)  at  Areclbo.  See  also  (Grebowsky  et.  al., 

1975) .  Studies  such  as  (Taylor  et.  al.,  1968)  have  shown 
that  nitrogen  Ions  are  also  a significant  constituent  at 
high  latitudes.  The  compositions,  however,  are  only  needed 
by  the  model  computer  program  inside  the  plasmasphere  (never 
at  high  latitudes)  so  that  ignoring  the  presence  of  nitrogen 
is  a justifiable  simplification. 

For  field  lines  outside  the  plasmasphere,  where  simple 
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diffusive  equilibrium  does  not  hold,  a convenient  model 
suggested  by  Angeram!  and  Carpenter  (1966)  Is  adopted.  This 
model  Is  expressed  by  (2.22)  below: 

M a Neo  (i.i*) 

Anger am  I and  Carpenter  (1966)  pointed  out  that  this  simple  r 
model  closely  approximates  the  more  accurate  Mcol 1 Is ionless" 
model  obtained  from  (Evlatar  et.  al.,  1964),  and  has  the 
advantage  of  extreme  simplicity. 

The  collisionless  model  assumed  a pure  hydrogen  plasma, 
while  observations  of  a light  Ion  trough  show  that  for 
ionospheric  heights  hydrogen  becomes  a minor  constituent  at 
high  latitudes  (Grebowsky  et.  al.,  1975).  The  error 
introduced  should  not  significantly  affect  the  conclusions 
of  this  analysis  because  the  stations  we  will  be  concerned 
with  are  seldom  poleward  of  the  plasmapause.  The  ray  paths 
more  frequently  encounter  the  collisionless  model  at  great 
heights,  beyond  the  plasmapause,  where  the  contribution  to 
n is  negligible. 

2. 1.3. 3 "Ionosphere" 

The  heart  of  the  semi empt ri cal  model  used  in  this  study 
is  N(h)  data  obtained  from  the  Isis  It  topside  sounder 
satellite  launched  April  1,  1971,  A sample  page  of  data  is 
shown  In  figure  11*2.  As  seen  from  the  sample,  the  soundings 
were  taken  about  every  fifteen  seconds,  or  about  every 
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Figure  11-2  A sample  page  of  Isis  II  Interpolated 


Ne(h)  data 
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degree  of  geographic  latitude.  The  quality  (4  In  figure 
11-2)  varies  considerably.  Best  quality  lonograms  are 
assigned  Q values  of  4 while  the  lonograms  with  the  least 
information  are  assigned  a Q value  of  9.  According  to  the 


introduction  to 

the 

compl 1 at Ion 

of  Isis 

1 1 data 

(see 

Acknowledgements ). 

for 

the  best 

qua  1 1 ty 

lonograms. 

the 

accuracy  in  the 

real 

height  is 

probably 

better  than 

5*. 

Unfortunately.  Ne(h)  Is  not  always  given  down  to  hmax. 
Furthermore,  no  Indication  in  the  data  compilation  Is  given 
where  this  occurs,  so  In  some  cases  hmax  is  actually  lower 
than  the  figure  we  are  forced  to  take. 

As  noted  earlier  (section  2. 1.3.1).  a latitude  vs. 
height  "real  density  array"  is  constructed  using  this  data. 
Linear  interpolation  is  used  between  the  data  points  to 
obtain  the  electron  densities  at  1000  and  700km  heights: 
NC1000)  and  NC700).  respectively.  The  other  parameters  read 
by  the  model  program  are:  Nmax.  hmax.  and  yzero.  where  yzero 
is  the  height  at  which  the  electron  density  of  the  model 
drops  to  zero  on  the  bottoms  I de.  Ne(h)  below  1000km  In 
height  is  treated  in  three  parts:  (a),  from  700km  to  1000km 

an  exponential  of  scale  height  HH  Is  matched  with  (b).  an 
exponential  of  scale  height  HL  between  hk  and  700km; 
finally,  (c).  a bi-parabola  of  semi thl ckness  ym  extends 
around  the  peak  from  y zero  to  hk.  where  hkahmax+ym/4.  and 
ymahmax-yzero.  The  equation  of  the  bi-parabola  is: 

h 


F 
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At  the  height  hk: 

N * Nma*  3 Nk  (i.A4) 

The  scale  height  HL  ts  obtained  from  (2.26)  by  rearranging 
the  exponential  equation  (2.25)  and  substituting 
HL,Nk,hk,N(700)  and  700  for  H,No,ho,N,  and  h respectively: 

N(W  * N.  «p (-hr*1-)  (»■«> 


HI  s 100  - hk tl.it) 

HL  IntNk/NOOO) 


Similarly  for  the  upper  exponential: 


131.317) 


[N(700)/tf(lW>0) 

Now  that  the  scale  heights  have  been  obtained,  the  computer 


program  can  calculate  N(h)  at  selected  heights.  A step  size 
of  20km  was  chosen  as  It  Is  small  enough  to  prevent  loss  of 
Important  detail  near  the  peak.  It  is  now  necessary  to 
consider  the  range  of  latitudes  over  which  the  ray  path 
passes  below  1000km.  This  range  is  divided  Into  51  vertical 
sections,  each  section  corresponding  to  one  of  the  51 
heights  along  the  ray  path  (again  the  step  size  is  20km) 
where  the  electron  density  is  calculated,  in  order  to  fill 
In  this  second  array  which  we  call  an  "interpolated  density 
array,"  1 1 wl 1 1 be  necessary  to  interpolate  between  the 
latitudes  of  the  real  density  array  and,  occasionally,  to 
extrapolate  beyond  the  range  of  the  real  density  array.  The 
Interpolation  Is  linear,  as  Is  the  extrapolation,  but  the 
number  of  points  used  to  determine  the  extrapolated  line 
depends  on  the  position  of  any  trough,  in  extrapolating 


i 


north  or  south  of  the  real  array  at  least  the  last  two 


points  are  taken.  Additional  points  are  taken  as  long  as 


each  additional  point  differs  from  the  preceding  one  by  less 


than  25$,  if  the  first  point,  on  the  other  hand,  differs  by 


at  least  25$  from  the  second,  a third  point  is  taken. 


Additional  points  are  taken  if  necessary  until  the  needed 


extrapolated  value  is  positive.  A least  squares  subroutine 


is  used  to  fit  a line  to  the  selected  points,  and  this  line 


provides  the  extrapolated  values  needed. 


Once  the  interpolated  density  array  has  been 


constructed,  the  vertical  TEC  at  51  columns  along  the  ray 


path  is  found  simply  from  the  vertical  segments.  The  slant 


TEC  Is  added  up  along  the  ray  path  by  multiplying  the 


electron  density  at  each  point  along  the  ray  path  by  the 


distance  between  points.  The  vertical  TEC  increment  along 


the  ray  path  is  then  obtained  by  multiplying  each  slant  TEC 


increment  by  cosX.  The  sum  of  these  vertical  increments  will 


be  called  here  the  “vertical  equivalent  TEC."  The  M factor 


is  then  calculated  at  each  point  along  the  ray  path  by  using 


(2.15)  and  the  Faraday  rotation  is  found  by  transforming 


(2.14)  into  a sum: 


xi  * 4fN(i)mi)Ah  l»# 

r i,i 


After  the  total  Faraday  rotation  along  the  ray  path  is 


computed,  it  is  combined  with  the  M value  at  each  height  to 


give  51  values  of  equivalent  vertical  TEC.  These  may  be 
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compared  with  the  TEC  of  the  vertical  columns. 

Tltheridge's  C1972)  conclusion  that  for  a wide  variety 
of  conditions  the  use  of  h(M)»420km  converts  the  Faraday 
rotation  to  TEC  up  to  2000km  leads  us  to  extend  our  vertical 
columns  up  to  2000km  so  that  we  may  compare  them  with  the  51 
equivalent  vertical  TEC  values  just  mentioned.  This  cannot 
be  done  to  great  accuracy  without  quadrupling  the  computer 
printout/  but  for  our  purposes  we  simply  integrate  from 
1000km  to  2000km  using  the  closed  expression: 

ffiodh  = HN(I000)[|  can) 

1 000  ^ 

A scale  height  of  H«1000km,  for  example/  simplifies  the 

reSU'ffiHA  = HN(iooo)(|--V)  - N(IOOO)  U3b ) 

II  9*0 

where  NC1000)  Is  In  units  of  10*e1/cm*  and  the  TEC  Is  In 
units  of  lO’^el/cm*-. 


2*2  Usefulness  of  constant  fl-factor  at  a high 
latitude  station  during  winter  night 

The  usefulness  of  equivalent  vertical  TEC  below  1000km 
except  in  the  presence  of  a steep  trough/  has  already  been 
shown  at  ml d-1 at! tudes  by  Haj eb-Hossel ni eh  (1973)*  We  extend 
his  work  to  high  latitudes  and  also  consider  the  region 
above  1000km  for  the  particular  case  of  winter  nighttime 
condl tions. 

2.2.1  Idealized  Arrays 

Consider  a uniform  data  array  extending  from  24*  to  67# 
geographic  latitude  with  the  following  Input  parameters: 
Nmax-2xl0*,  N(700)-2. 63x10*,  N(l000)-lxl04,  hmax-320km,  and 
yzero«200km.  The  ray  path  travels  from  a geostationary 
satellite  to  each  of  three  stations,  the  coordinates  of 
which  are  found  in  Table  11*1. 


Table 

1 l-l 

station 

geog.  iat. 

geog.  long. 

m-'g.  1 

Haml 1 ton 

42.63 

70.82W 

54.34 

Goose  Bay 

53.30 

60.33W 

64.82 

Narssar ssuaq 

61.20 

45.40W 

71.38 

Consider  the  input  data  shown  in  Table  11*11  (the 
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Table  11-11 


Uniform  Latitude  Gradient  Input  Array 


geog, 

lat. 

N(1000)  „ 
(10* 

NC700) 
el/cnr  ) 

Nmax 

hmax 
(km. ) 

yzero 
(km. ) 

25 

1.6 

3.5 

14.0 

320 

200 

28 

1.5 

3.3 

13.4 

32  0 

200 

31 

1.4 

3.1 

12.8 

320 

200 

34 

1.3 

2.9 

12.2 

320 

200 

37 

1.2 

2.7 

11.6 

320 

200 

40 

1.1 

2.5 

11.0 

320 

200 

43 

1.0 

2.3 

10.4 

320 

200 

46 

0.9 

2.1 

9.8 

320 

200 

49 

0.8 

1.9 

9.2 

320 

200 

52 

0.  7 

1.  7 

8.6 

320 

200 

55 

0.6 

1.5 

8.0 

320 

200 

58 

0.  5 

1.3 

7.4 

520 

200 

61 

0.4 

1.1 

6.8 

320 

20  0 

64 

0.3 

0.9 

6.2 

320 

200 

67 

0.2 

0.7 

5.6 

32  0 

200 

70 

0.1 

0.5 

5.0 

320 

200 

densities  are  In  units  of  loVcm*  ). 

It  represents  a uniform  latitude  gradient.  Examination 
of  Brace  et.  ai.  C1967)  will  show  that  the  data  was  chosen 
as  a reasonable  compromise  between  average  winter  nighttime 
conditions  and  the  requi renents  of  constructing  a uniform 
density  gradient  over  a large  range  of  latitudes. 

Examination  of  the  "vertical  equivalent"  TEC 
accumulated  for  each  of  these  two  arrays  exclusively  using 
the  diffusive  equilibrium  model  for  heights  above  1000km 
leads  to  the  the  conclusion  that  a scale  height  H«1000km  in 
(2.29)  will  reproduce  the  results  to  the  required  accuracy. 
Using  (2.30)  we  now  extend  the  51  vertical  columns  to  2000km 
and  compare  with  the  equivalent  vertical  TEC  computed  for 
the  51  columns.  The  results  are  shown  in  Table  11*111  where 
"best  h(R)"  refers  to  the  height  at  which  the  vertical 
columnar  TEC  matches  the  equivalent  vertical  TEC  using  the  R 
factor  for  that  height.  "%dlf  at  42d"  gives  the  percentage 
difference  between  the  vertical  TEC  at  the  420-km  point  and 
the  equivalent  vertical  TEC  using  the  R factor  at  420km. 

Table  11*111 


station 

best 

h(R) 

?dl  f 

at  420 

Uni  f 

Uni  f 

Uni  f 

Uni  f 

data 

grad 

data 

grad 

Hami 1 ton 

371km 

220  km 

2.8 

7.4 

Goose  Bay 

400km 

1 7km 

1.2 

9.1 

Narssar ssuaq 

33  3km 

4 km 

2.1 

17.7 
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The  results  clearly  show  that  the  use  of  h(R)=420km 
gives  good  agreement  for  a uniform  data  array.  For  a 
gradient/  however/  there  is  decreasing  agreement  with 
increasing  station  latitude.  Caution  must  be  used  in 
considering  the  results  since  both  of  these  arrays  are  quite 
i deal ized. 

It  is  the  gradient,  and  not  a uniform  data  array  which 
more  closely  approximates  typical  winter  nighttime 
conditions.  Examination  of  the  computer  output  shows  that 
the  Faraday  rotation  technique  as  far  poleward  as  the  trough 
minimum  apparently  measures  total  electron  content  not  to 
2000km/  but  to  some  greater  height  for  these  winter 
nighttime  conditions.  The  effect  is  due  to  the  fact  that 
heights  along  the  ray  path  in  the  range  2000km  to  5000km 
coincide  with  field  lines  of  the  lowest  invariant  latitude/ 
i . e.  / the  the  highest  reference  level  (1000km)  electron 
densities.  This  enhances  the  Faraday  rotation  accumulated  in 
this  upper  region.  Furthermore,  the  minimum  invariant 
latitude  corresponds  to  higher  heights  for  higher  latitude 
stations,  enhancing  the  effect  there.  For  the  uniform 
gradient  the  minimum  invariant  latitude  occurs  at  about 
2500km  observing  from  Hamilton,  Mass.  (4  2°.  ON,  70°.  8W  ),  3700km 
from  Goose  Bay,  Labrador  ( 53°.  3N,  bu°.  3W),  and  4200km  from 
Nar  ssarssuaq,  Greenland  (61°,  2N,  45*.  **W).  The  geometry  is  shown 


In  figure  11-3,  the  lowest  invariant  latitude  the  ray 
path  crosses  at  one  point  (point  B),  corresponding  to  the 
height  h 1.  Points  A and  C have  the  same  invariant  latitude 
and  correspond  to  smaller  and  greater  heights  respectively. 
For  a lower  latitude  station  the  lowest  invariant  latitude 
corresponds  to  a lower  height  h2‘. 

To  check  that  this  effect  is  real  and  not  due  to  the 
use  of  an  idealized  gradient,  a real  gradient  taken  from 
Isis  II  topside  data  was  used  (December  15,  ld71).  The 


conclusions  above  were  confirmed. 

Clearly  a non-negl i gi bl e error  is  introduced  in 
assuming  that  the  W value  used  to  calculate  equivalent 
vertical  TEC  up  to  lOUJkm  is  little  different  from  the  ffi 
value  used  to  calculate  equivalent  vertical  TEC  up  to  20i)0km 
or  higher.  For  the  uniform  gradient  observed  from  Goose  Bay, 
for  instance,  17.5*  of  the  "vertical  equivalent  TEC"  is 
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accumulated  between  1000km  and  3700km,  the  upper  height 
being  the  height  down  to  which  the  first  10®  of  Faraday 
rotation  (a  measurable  amount)  was  accumulated.  The  R factor 
used  to  compute  the  TEC  up  to  1000km  will  then  differ  by 
17. 5i  from  the  ffl  factor  used  to  compute  the  TEC  up  to 
3700km. 

In  figure  11-4,  the  hatched  regions  refer  to  distances 
along  ray  paths  to  Hamilton  and  Goose  Bay  where  60%  of  the 
total  Faraday  rotation  (0*,*)  *s  accumulated,  with  10%  of 
above  and  10%  of  DToT  below. 

2.2.2  Severe  Latitude  Gradient 

Consider  the  data  array  shown  in  Table  1 1 — I V . 

This  data  was  obtained  between  2:19  UT  and  2:30  UT  on 
December  18,  1971  when  a severe  magnetic  storm  was 

occurring.  At  the  time  the  data  was  taken  the  Kp  index  stood 
at  4-.  Six  hours  earlier  it  had  reached  7+.  The  parameters 
Nmax,  N ( 7 0 0 ) , N(1000),  and  N(sat)  are  plotted  in  figure 
11-5.  N(sat)  Is  the  measurement  closest  to  the  satellite, 
always  near  1400km.  One  can  clearly  see  the  sharp 

equatorward  edge  of  the  trough  near  40®  which  is  identified 

» i 

with  the  plasmapause.  The  poleward  edge  of  the  trough  is 
also  steep  and  north  of  the  trough  the  profile  is  erratic. 

The  yzeros  were  somewhat  arbitrarily  assigned  but  the 
presence  of  the  auroral  E layer  near  the  trough  poleward 

lL _ 
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Table  I I -IV 

Severe  Latitude  Gradient  Input  Array 


December  18,  1971 


geog. 
1 at. 

geog. 

long. 

N(1000) 

(10* 

N( 700 ) 
el  / cm5 

Nmax 

) 

hmax 
( km,  ) 

yz  ero 
(km. ) 

30.09 

53.  10 

1.37 

4.45 

16.26 

520 

195 

3d . 16 

52.80 

1.13 

3.7  S 

12.51 

500 

195 

39.  90 

52.  80 

. 5d 

1.24 

3.94 

300 

195 

42.12 

52.  70 

1.  00 

2 , 2u 

8.80 

320 

195 

43.65 

53.  70 

1.05 

2.65 

8.77 

3*+0 

195 

46.67 

52.  70 

1.34 

3.  76 

11.04 

3u0 

195 

48.  15 

52.  70 

1.  53 

4.24 

9.92 

420 

195 

48.  88 

52.60 

2.81 

7.85 

17.51 

*480 

165 

49.62 

52.60 

1. 67 

4.8  7 

12.  23 

400 

165 

50.  41 

52.60 

2.50 

o . 0 1 

19.00 

340 

lo  5 

51.21 

52.60 

. 97 

2.32 

12.23 

260 

163 

52.68 

52.  70 

.48 

1.11 

15. 6o 

200 

120 

54.  21 

52.  70 

.32 

.80 

0.57 

20U 

120 

55.68 

52.  70 

. 31 

. 80 

16.01 

20  0 

120 

60.  16 

52.80 

. 32 

. 78 

7.90 

260 

165 

67.  35 

53.50 

. 15 

. 60 

2.35 

360 

165 
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edge  necessitates  the  low  yzero  values  there. 

To  see  the  effect  of  this  severe  trough  on  the 
calculation  of  TEC/  we  have  the  ability  to  shift  the  data 
array  by  any  number  of  degrees  north  or  south  observing  it 
from  any  given  station.  Goose  Bay  is  chosen  (near  L»4)  and 
the  array  is  shifted  north  and  south  from  1U • to  -5#  in  one 
degree  steps.  In  Table  ll-V  the  results  are  shown.  The 
height  or  heights  h CR ) along  the  ray  path  at  which  the 
vertical  columnar  TEC  matches  the  equivalent  vertical  TEC 
obtained  from  the  M value  at  that  height  or  heights,  and  the 
percentage  difference  between  vertical  TEC  at  the  420km 
point  and  the  equivalent  vertical  TEC  using  the  R value  at 
42Ukm  are  listed  for  each  shift. 

For  a given  shift  there  may  be  up  to  four  heights  where 
equivalent  vertical  TEC  and  vertical  TEC  match.  This  renders 
selection  of  a best  height  h(M)  Impossible.  The  average 
difference  between  vertical  and  equivalent  vertical  TEC  at 
the  420km  point  is  26%,  much  larger  than  the  5%  variation 
usually  tolerated.  One  also  must  expect  very  large 
variations  on  occasion  with  such  a steep  trough  as  indicated 
by  the  8°  north  shift's  151%  variation.  To  see  why  such  a 
large  discrepency  appeared  for  the  + 8*  shift  the  electron 
density  along  the  ray  path  below  1000km  has  been  plotted  in 
figure  1 1 -6.  Because  the  ray  path  passes  through  the  wall  of 
the  trough  equatorward  edge  well  above  hmax  of  the  vertical 
column  the  profile  is  greatly  distorted.  Since  the  vertical 
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Table  I l-V 

. 


shift  best  h(H)  % dif  at  420km. 


♦ 10 

none 

6.5 

♦9 

9 

4.8 

♦ 8 

67,  430 

151. 

♦7 

940,  557,  136 

73. 

♦6 

331,  641 

15. 

♦ 5 

397, 740 

2.2 

♦4 

426, 850 

1.1 

♦3 

462 

6.7 

♦2 

508 

13. 

♦1 

539 

15. 

0 

161,  403 

1.0 

-1 

267, 360, 456 

25. 

-2 

359,  412,497,  527 

5.1 

-3 

429, 515, 548,630 

9.3 

-4 

504,  620,  73a 

41. 

-5 

596,  720,834 

49. 

column  is  inside  the  trough  the  vertical  total  content  is 
far  less  than  the  "vertical  equivalent"  TEC. 

It  appears  that  for  a severe  trough  one  cannot 
legitimately  convert  Faraday  rotation  into  the  total 
electron  content  of  a vertical  column/  confirming 
Haj eb-Hossei nl eh  (1973). 

2.2.3  Concl us i ons 

For  the  winter  nighttime  conditions  we  have  examined  at 
high  latitudes,  Faraday  rotation  apparently  is  a measurement 
of  equivalent  vertical  TEC  up  to  some  height  higher  than 
2000km.  A more  refined  model  is  needed  to  fix  the  actual 
height.  Under  conditions  of  a severe  trough  Faraday  rotation 
cannot  be  meaningfully  converted  to  total  electron  content 
using  any  mean  h(R).  In  attempting  to  understand  such 
conditions  it  Is  best  to  study  the  variations  in  Faraday 
rotation  without  attempting  to  convert  it  into  the  total 
content  of  a vertical  column.  Further  study  is  needed  to 
determine  the  dependence  of  the  best  h(M)  on  seasonal,  solar 
cycle,  diurnal  variations,  and  differences  In  station 
latitude  in  the  absence  of  a severe  trough. 


2.3  lonosondes  and  Slab  Thickness 


2.3.1  lonosondes 


The  Isis  II  satellite  observations  described  earlier 
were  taken  by  an  ionosonde,  lonosondes  collect  huge 
quantities  of  data  from  the  ground  as  well/  and  in  fact  are 
the  most  widely  used  means  of  monitoring  ionospheric 
electron  density  in  addition  to  many  other  ionospheric 
parameters.  Ionosonde  observations  taken  from  the  surface 
of  the  Earth  probe  the  "bottoms I de"  ionosphere/  l.e.,  the 
ionosphere  below  hmax,  where  the  last  reflection  occurs. 
Similarly/  satellite  sounders  such  as  that  on  Isis  II  probe 
the  "topside"  ionosphere. 

The  ionosonde  makes  use  of  pulsed  radio  signals/  each 
pulse  sweeping  over  a range  of  frequencies  from  1-2  MHz  to 
perhaps  25  MHz.  The  receiver  displays  on  an  oscilloscope  a 
trace  showing  the  frequency  of  the  returned  signal  vs.  time. 
The  virtual  height  h'  at  which  the  signal  was  reflected  is 
then  gi ven  by: 


I 


h*  - ict  u.*i) 


In  order  to  see  why  the  virtual  height  and  the  true 
height  of  the  reflection  differ  it  is  necessary  to  review 
the  process  of  reflection.  For  this  simple  argument  we 
Ignore  the  presence  of  the  magnetic  field  and  of  collisions. 
Then  the  transmitted  radio  wave  passes  through  a medium  with 
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an  Index  of  refraction  given  by  (2.6)  with  Y^«0,  l.e,. 


■FW 


The  phase  velocity  of  the  wave  Is  then  related  to  the  index 
of  refraction  and  the  group  velocity  by  the  next  two 
relations : 


y * cl y (2.33) 

\^Vy»  c*  (2.34) 

As  the  radio  wave  goes  Into  regions  of  higher  electron 
density  the  plasma  frequency  increases  by  (1.1)  which 
decreases  yj,  Increases  v^  and  decreases  Vy  thus  slowing  the 
wave  packet  down.  When  * f the  group  velocity  becomes 
zero  and  the  wave  Is  reflected.  Thus  we  see  that  the 
reflection  Is  not  instantaneous  but  rather  gradual  so  that 
the  virtual  height  Is  always  greater  than  the  true  height  of 
the  reflection.  The  relation  between  the  two  is  given  by  the 
process  of  true  height  analysis,  too  complicated  to  be  even 
outlined  here.  For  a brief  description  see  (Rishbeth  and 
Garrlott,  1969). 

True  height  analysis  is  costly  and  not  often  necessary. 
The  peak  density  NmF2  can  be  easily  obtained  from  (1.1)  by 
reading  foF2  right  off  the  ionogram  trace.  foF2  is  always 
obtained  from  a vertical  reflection.  The  maximum  frequency 
reflected  to  a 3000km  distance  is  known  as  MUF(3000)F2.  An 
empirical  formula  known  as  the  Shimazaki  relation  gives  a 
good  approximation  to  the  true  height  of  the  reflection, 
especially  at  night  (Shimazaki,  1955)  using  MUF(3000)F2  (or 
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M(3000)F2).  It  is: 

l c*\  - - 176km 

nmr<*  - M(3ooo)F3. 


(2.3$) 


While  refinements  have  been  proposed  by  Wright  and  McDuffie 
(1960)  and  by  Bradley  and  Dudeney  (1973),  (2.35)  remains 
fairly  accurate  for  nighttime  conditions. 


2,3.2  Slab  Th i ckness 


Slab  thickness  is  defined  as  the  ratio: 


T = TEC/Nmax 


(2.36) 


and  thus  represents  the  thickness  the  ionosphere  would  have 
if  the  electron  density  throughout  was  the  constant  peak 
value.  Here  Nmax  must  be  the  peak  density  of  the  vertical 
TEC  column  unless  there  are  very  shallow  horizontal 
gradients  In  TEC  and  Nmax. 

Slab  thickness  is  an  Important  ionospheric  parameter 
since  monitoring  it  gives  an  indication  of  the  shape  of  the 
vertical  profile  of  N(h)  without  the  need  to  resort  to  true 
height  analysis.  If  the  total  content  of  a vertical  column 
remains  constant  while  the  peak  density  changes,  for 
example,  one  must  conclude  that  a redistribution  of  plasma 
has  taken  place,  i.e.,  the  shape  of  the  profile  of  electron 
density  vs.  height  has  changed.  Alternately,  if  TEC  changes 
and  Nmax  does  not,  one  concludes  that  an  influx  of  plasma 
has  occurred  probably  bringing  electrons  down  from  the 
exosphere  into  the  vertical  column.  As  seen  In  section  2.2 
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one  must  be  cautious  about  making  such  conclusions  using 
Faraday  rotation  measurements  in  the  presence  of  troughs 
since  the  equivalent  TEC  and  not  the  actual  TEC  of  the 
column  is  known. 

In  chapter  four  we  will  see  examples  in  which  slab 
thickness  decreases  significantly  In  the  presence  of  troughs 
on  the  basis  of  Faraday  rotation  measurements.  We  will 
Investigate  the  possible  roles  of  both  trough  motions  and 
vertical  distortions  in  producing  this  phenomenon. 
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3.1  Mean  and  Normalized  Monthly  Curves 
Median  Behavior 

In  the  previous  chapter  we  discussed  the  importance  of 
three  ionospheric  parameters:  total  electron  content  (TEC), 
peak  density  of  the  F-2  layer  (Nmax),  and  the  ratio  of  these 
two  quantities,  the  equivalent  slab  thickness  (#t*TEC/Nmax) . 

The  purpose  of  this  chapter  is  to  briefly  discuss  the 
average  seasonal  and  diurnal  behavior  of  these  three 
parameters  at  60  degrees  magnetic  latitude. 

In  this  chapter,  the  term  "TEC,"  will  refer  to  the 
equivalent  vertical  TEC  through  the  420-km  point,  i.e.,  in 
order  to  obtain  the  TEC,  Faraday  rotation  measurements  are 
combined  with  a constant  R value  at  420km  height  along  the 
ray  path.  With  TEC  defined  In  this  manner,  the  terms  total 
electron  content  (TEC)  and  Faraday  rotation  (ft)  may  be  used 
i nterchangeabl y. 

38  months  of  TEC  data  collected  at  Goose  Bay,  Labrador 
( 53® 3geog. N, 60? 3geog. W)  will  be  combined  to  calculate  the 
mean  monthly  behavior  of  TEC.  A vertical  column  intersecting 

I 

the  ray  path  to  the  geostationary  satellite  (ATS-3)  at  420km 
height  is  located  at  geographic  coordinates  47?5N,  62*.  2W,  and 
generalized  invariant  latitude  59?9  (L-3.98).  Fortunately, 
we  have  available  many  years  of  N max  data  taken  by  the 
tonosonde  at  St.  John's,  Newfoundland  (47%geog.N, 
52?7geog.W),  located  at  nearly  the  identical  latitude  as  the 
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420-km  point  and  only  thirty  minutes  of  local  time  to  the 
east.  41  months  of  St.  John's  Nmax  data  spanning  the  same 
period  as  the  TEC  measur&nents  will  give  us  the  mean  monthly 
Nmax  behavior.  Then  we  will  combine  the  TEC  and  Nmax  data 
into  mean  monthly  slab  thickness  results  based  on  57  months. 
The  months  from  which  the  mean  monthly  behavior  of  the  three 
parameters  are  derived  are  listed  in  Table  lll-l. 

Table  lll-l 

Data  Base  for  Mean  Monthly  Calculations 

Parameter 


TEC 

Nov. , 

l'J  7 1-No v. , 19  73;  Apr  l 1 , 

1974-Apr  i 1, 

1975 

Nmax 

Nov. , 

1971-March,  1975 

* 

Nov. , 

1971-Nov.,  1973;Apr  i 1, 

19  74-March, 

1975 

The  monthly  mean  diurnal  behavior,  which  for  TEC  we 
call  MTEC(  Y,  M,  H ),  Is  found  by  first  calculating  the  mean 
values  at  each  hour  of  local  time  for  the  3 parameters 
during  each  month.  This  gives  24  numbers  for  each  month.  All 
the  Januarys,  Februarys,  etc.,  are  then  combined  by  taking 
the  mean  of  the  mean  values  for  each  local  time.  This  gives 
mean  monthly  diurnal  curves  which  for  TEC  we  call 
MMTEC(M,  H),  and  also  curves  for  Nmax  and  *t.  These  curves  are 
presented  in  figures  lll-l,  I I 1-2,  and  II 1-3. 

In  calculating  these  mean  monthly  diurnal  curves,  we 
have  followed  (Hawkins  and  Klobuchar,  19  74).  In  that  study, 
seven  years  of  TEC  data  collected  at  Hamilton,  Mass.  (42*6 
geog. N, 7Q?8geog. W)  were  presented.  We  will  be  interested  in 
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Figure  MI-3  Mean  monthly  slab  thickness  diurnal  curves 
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comparing  the  average  behavior  found  for  the  mid-latitude 
ionosphere  (L»2.8)  in  that  study  with  the  l«4  data  presented 
here  (see  sec.  3.2). 

Hawkins  and  Klobuchar  (1974)  noticed  that  the  highest 
noontime  TEC  values  at  ini  d 1 at  I tudes  occurred  near  the 
equinoxes  (March  and  October).  This  effect  is  attributed  to 
relative  changes  in  the  concentrations  of  atomic  oxygen  and 
molecular  nitrogen  (the  0/N4  ratio),  which  are  responsible 
for  the  producti on/ loss  process  in  the  F-region.  We  notice 
that  the  highest  mean  monthly  TEC  at  1=4  occurs  during  the 
daytime  of  winter  and  equinox  months  (figure  lll-t).  In 
general,  however,  these  "seasonal"  and  "equinoctal" 
anomalies  do  not  seem  to  be  as  prominent  in  the  mean  Goose 
Bay  results  as  In  those  from  Hamilton.  Of  course,  the  Goose 
Bay  data  refer  to  higher  latitudes  and  the  data  base 
represents  years  of  low  solar  flux.  It  appears  that  under 
these  conditions,  i.e.,  L=4  and  low  sunspot  number,  such 
anomalies  in  the  mean  diurnal  curves  tend  to  decrease  in 
promi nence. 

In  order  to  eliminate  this  effect,  and  also  to 
eliminate  the  effect  of  variations  in  solar  flux,  HawKins 
and  Klobuchar  (1974)  normalised  the  monthly  TEC  values.  We 
follow  their  procedure. 
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First  the  monthly  mean  TEC  values  (MMTEC)  for  a period 
of  24  hours  are  averaged  together: 

ATEC(Y,M)  = pi  ^ MTECCY.M.H)  (3.1) 

Then  normalized  monthly  mean  curves  are  obtained  via: 
NTEC(Y,M,H)  * MTEC(Y,M,H)/ATEC(Y,M)  (3.2) 

The  normalized  curves  were  averaged  over  the  years  of  the 
study  by  calculating: 

, Y' 

NMTEC(M.H)  = -Jr  £ NTECAy.M.H) 

Y«l 


These 

mean 

normalized  monthly  TEC 

curves 

remove  any 

var i at i on 

due 

to  the 

solar  flux  and 

the 

"equi noctal 

anomaly,  11 

1 arge 

effects 

from  wh i ch  wou Id 

tend 

to  increase 

ATEC(Y, M) . 

What 

r emai ns 

is  the  seasonal 

variation  in  the 

diurnal  curves.  The  results  of  this  analysis  for  TEC  appear 
in  figure  111-4.  The  same  procedure  has  been  followed  for 
Umax  and  't,  and  the  normalized  curves  for  these  parameters 
are  found  in  figures  111-5,  and  1 1 I — o . 

Note  the  smooth  progression  in  noon-time  normalized  TEC 
(figure  111-4),  with  highest  values  in  winter,  and  lowest  in 
summer.  These  curves  show  that  the  winter  ionosphere  is 


solar  production  dominated,  responding  primarily  to  changes 
in  solar  zenith  angle.  The  seasonal  anomaly  is  strikingly 
shown,  with  the  lower  summer  daytime  values  due  to  enhanced 


i 


Figure  111-6  Normalized  mean  monthly  slab  thickness 


diurnal  curves 
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plotted  in  figures  II 1-7  and  111*8,  respectively. 

In  figure  111*8,  curves  for  both  Hamilton  (from  the 
Hawkins  and  Klobuchar  study)  and  Goose  Bay  are  provided. 
Notice  how  the  smooth  seasonal  progression  apparent  for  the 
maximum  values  in  figure  111-8  is  not  found  in  figure  III-/. 
The  minimum  values  follow  a fairly  regular  progression  in 
both  figures.  One  thing  of  interest  is  the  close  correlation 
between  the  curves  for  Hamilton  and  Goose  Bay  in  figure 
111*8. 

Monthly  conditions  for  the  TEC  observed  from  Goose  Bay 
are  presented  in  the  form  of  a contour  plot  in  figure  II 1-9. 
This  plot  covers  all  38  months  of  TEC  data.  Times  of  sunrise 
and  sunset  at  the  420-km  point  are  indicated  oy  tics. 

We  will  discuss  the  general  features  of  ionospheric 
behavior  at  L»4  In  the  next  section  with  reference  to  the 
curves  already  presented. 
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MAXIMUM  AND  MINIMUM 

MEAN  MONTHLY  TEC  . 

GOOSE  BAY  (1071-1075) 

MONTH  J/IN  FEB  MAR  APR  W OEM  -JUL  AUG  SEP  CCT  UOV  DEC 

DAY  NO.  IS  45  74  105  IBS  ICC  1%  117  S 123  3R  34A 

HIGHEST  IC.3  I8.C  17.5  15.4  173  lfc.4  15.1  14.7  IC.5  I P.7  ICS  17.1 

LOWEST  l.b  13  3.1  3.3  4.4  4.C  4.0  3.0  11  14  1.1  1.3 


MAXIMUM  AND  MINIMUM  MEAN  MONTHLY 

NORMALS  ED  TEC 
HAMILTON  UU7-IW 

MONTH  m FEB  MAR  APR  MAY  JUN  JUl  UK,  SEP  OCT  NCV  DEC 

DAY  NO  15  45  74  105  135  155  1%  527.13  188  318  348 

HIGHEST  .5.10  1.45  1.77  1.55  1-50  (.45.  (.44  1.48  1.51  HI  1.05  110 

LOWEST  JS  ..30  .53  . 35  . 33  .41  .38  ,3a  .IT  .15  .31  .33 


GOOSE  BAY  (H71- 

HI6HEST  _3.331.I5  1.84  1.53  1.58  1.50  |.51  1.54  1.30  1.97  1.18  1.37 
LOWEST  .55.  .11  .31  .34  .33  .41  .40 ...33  .34  .55  .18  .35 


MAX  IMUL 


V'l  ■ ; DAY  NUMBER 

* FfiPMCHaw^inS  anJ  iClotuckar,  H74) 


GOOSE  BAY  MEDIAN  TEC  (10**1  2 EL/  CM**2  ) 
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3.2  Discussion 

We  now  wish  to  briefly  discuss  the  general 
characteristics  of  the  ionosphere  as  observed  from  Goose  Bay 
(TEC  measurements)  and  St.  John's  (Umax  results)/  and  to 
compare  these  features  with  available  midlatitude  results. 

Figure  I I I — I has  shown  that  the  general  features  of  the 

mean  monthly  TEC  at  Goose  Bay  are  the  sane  as  those  at 

Hamilton,  which  Hawkins  and  Klobuchar  (1974)  described  as 
fol lows : 

"a.  A decay  during  the  hours  of  darkness, 

prolonged  and  1 i near  in  summer  months,  more 

abrupt  in  winter  with  a low-level  steady 
component  in  the  latter  part  of  the  night. 

"b.  An  abrupt  rise  at  dawn  at  approximately 
U900UT,  and 

"c.  A daytime  maximum  in  the  p.m.  hours 
with  a growth  and  decay  that  approximates  to 
a sinusoid  in  winter  and  a saw-tooth  in  summer." 

Comparison  of  the  mean  curves  at  Goose  Bay  and  Hamilton 
(figure  I I 1-1  here  and  figure  2 in  Hawkins  and  Klobuchar's 
(19/4)  study),  shows  that  for  each  month  and  each  hour,  mean 
monthly  TEC  observed  at  Goose  Bay  is  less  than  that  observed 
at  Hamilton.  This  is  to  be  expected  from  the  greater  solar 
zenith  angle  at  all  hours  at  the  Goose  Bay  420-km  point 
compared  to  the  Hamilton  42d-km  point,  and  from  the  higher 
sunspot  number  during  the  midlatitude  study. 

One  also  notices  that  the  seasonal  variation  in  maximum 
TEC  seems  to  be  less  at  Goose  Bay  than  at  Hamilton.  Since 
the  normalized  curves  (see  especially  figure  111-8)  show  a 
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nearly  identical  seasonal  variation  this  effect  must  be  due 
to  the  factors  mentioned  in  the  last  section#  viz.#  higher 
latitude  and  low  sunspot  number. 

The  evening  anomaly#  i .e.#  maximum  values  of  TEC  In 
summer  months  occurring  during  the  early  evening  hours#  Is 
certainly  present  both  In  the  Goose  Bay  and  in  the  Hamilton 
data.  In  the  Goose  Bay  results#  maximum  daily  TEC  occurred 
from  May  through  August  at  1800LT. 

Since  the  (Hawkins  and  Klobuchar#  1974)  study  did  not 
present  Nmax  and  t*  curves#  we  are  unable  to  compare  our 
results  for  St.  John's  Nmax  with  comparable  results  for 
Wallop's  Island  (or  Ft,  Belvoir)#  near  the  Hami i ton  420-km 
point.  We  can  notice#  however#  that  the  mean  monthly  Nmax 
curves  seem  to  be  similar  in  over-all  shape  to  the  mean 
monthly  TEC  curves.  The  seasonal  anomaly  (see  sec.  1.1.4)  is 
present#  as  Is  an  evening  anomaly.  For  the  evening  anomaly# 
we  notice  that  maximum  daily  Nmax  occurs  at  1900LT  from  May 
through  August. 

While  mean  monthly  mean  curves#  t.e.#  curves 
representing  several  years  of  data#  are  not  available  for 
mid-lati tudes#  Klobuchar  and  Allen  (1970)#  presented  monthly 
mean  slab  thickness  diurnal  curves  for  several  months  in 
1968.  From  these  curves  it  appears  that#  at  midlatitudes#  *t 
Is  lower  at  night  except  in  winter.  These  curves  also  show 
that  the  highest  noontime  't  values  occur  In  the  summer 
months#  with  a steady  progression  from  high  values  in  summer 
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to  low  values  in  winter.  This  latter  variation  is  also 
apparent  in  the  mean  monthly  ir  curves  for  Goose  Bay  and  St. 
John's  (figure  ill-3). 

As  in  the  case  of  midlatitudes,  there  is  no  clear 
monthly  progression  at  night,  although  nighttime  seems  to 
be  higher  in  winter  than  in  summer.  It  appears  that  *t  is 
higher  at  night  for  most  of  the  year,  contrary  to  the 
midlatitude  results. 

These  results  are  in  accord  with  our  generalization 
that  the  electron  density  trough  is  more  prominent  in  winter 
and  at  higher  latitudes,  because  while  Nmax  will  decrease  to 
very  low  values  in  the  trough,  the  Faraday  rotation,  which 
is  accumulated  along  a long  ray  path  passing  through  higher 
density  regions  to  the  south,  will  not  decrease  as  much. 
Therefore,  t^TEC/Nmax  will  increase  under  trough  conditions. 
Since  electron  densities  Increase  along  the  ray  path  south 
of  the  vertical  column,  the  Faraday  rotation  measures  the 
total  electron  content  of  the  vertical  column  to  heights 
greater  than  2000km  during  trough  conditions.  This  is  the 
same  conclusion  we  reached  in  Chapter  II  through  analysis  of 
a model  density  gradient. 

Let  us  now  turn  our  attention  to  the  normalized  curves 
(figures  1 1 1 — 4 to  ill-6).  Figure  ill-4  for  TEC  should  be 
compared  with  figure  8 in  (Hawkins  and  KloOuchar,  iy74). 
These  sets  of  curves  are  remarkably  similar.  The  seasonal 
progression  of  the  curves  is  Identical,  showing  a strong 
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seasonal  anomaly  for  both  sets.  The  close  correlation 
between  these  two  sets  of  curves  has  already  been  pointed 
out  by  reference  to  figure  I 1 1-8,  where  the  only  significant 
difference  seems  to  be  a slightly  greater  maximum  normalized 
TEC  for  all  seasons  at  L®4  than  at  mi dlati tudes.  This  means 
there  Is  a greater  difference  between  daytime  and  nighttime 
TEC,  which  is  apparently  due  to  very  low  values  at  night 
associated  with  the  electron  density  trough  in  the  Goose  Bay 
data. 

The  normalized  Nmax  curves  (figure  lll-a)  show  a very 
striking  seasonal  anomaly,  more  prominent  than  that  found 
for  the  mean  monthly  Nmax  (figure  1 1 1-2). 

The  normalized  slab  thickness  curves  (figure  I I 1-6) 
show  less  variation  from  month  to  month  than  the  mean 
monthly  X curves,  but  otherwise  show  basically  the  sane 
seasonal  progression,  with  lowest  daytime  X in  the  winter 
and  greatest  nighttime  X in  that  season. 

In  the  contour  plot  of  monthly  median  TEC  (figure 
I I 1-9),  which  covers  all  38  months  of  data,  we  can  clearly 
note  several  of  the  effects  previously  mentioned.  An 
"equinoctal  anomaly"  is  present,  as  seen  by  the  high  daytime 
values  in  March  and  October,  1972,  October,  1973,  and 
October,  1974.  There  is  unquestionably  a strong  variation 
f row  year  to  year.  Compare,  for  example,  March,  1972,  1973, 
♦ '•u  1*75.  The  maximum  daytime  value  in  March,  1972  is  over 
» . . „ *•* 1 /cm*.  In  March,  1973,  it  falls  to  about  17xl014, 
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and  In  March/  1975/  it  has  reached  as  low  as  9xlOia'.  The 
unfortunate  four  month  gap  in  the  data  prevents  us  from 
comparing  a March/  1974  value  with  the  values  for  the  other 
three  years.  This  yearly  variation  is  clearly  due  to  changes 
in  solar  flux.  The  mean  value  of  the  observed  solar  flux  at 
2800MHz  prepared  by  Ottawa  is  128.5  in  March,  19  72,  but  only 
72.4  in  March,  19/5.  As  1975  is  near  the  low  point  of  the 
solar  cycle,  the  rate  of  production  of  electrons  is  very  low 
at  that  time.  Other  things  to  notice  are  the  higher 
nighttime  values  of  TEC  In  summer,  and  the  presence  of  the 
evening  anomaly  in  that  season  (maximum  values  in  the  early 
eveni ng) . 
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3.3  Conclusions 

The  ionosphere  at  L=4  shows  many  similarities  with  the 
midlatitude  ionosphere.  Comparison  of  Goose  day  TEC 
(refering  to  L*4)  for  the  period  1J71-1J75  wi th  Hamilton  TEC 
(refering  to  L®2.tf)  for  the  period  lyo7-XB73  revealed  that 
the  Goose  Bay  data  showed  lower  TEC  on  the  average  for  al 1 
seasons  and  times  of  day.  This  is  due  to  two  factors:  the 
larger  solar  zenith  angle  at  all  times  at  L=4  in  comparison 
with  L=2.8,  and  the  fact  that  the  years  of  the  Hamilton 
study  were  years  of  higher  sunspot  number  than  the  years  of 
the  Goose  Bay  data.  The  other  differences  between  the  two 
data  sets  are  not  as  striking  and  can  all,  it  seems,  be 
traced  to  the  presence  of  the  electron  density  trough  at 
L=4,  most  prominent  during  winter  nights. 


i 


mi 


74 


4.1  Nighttime  Slab  Thickness  Decreases  Observed 
from  Goose  Bay  during  December,  1971 


As  discussed  in  the  previous  chapter,  total  electron 
content  obtained  at  Goose  Bay,  Labrador  may  be  combined  with 
simultaneous  peak  density  observations  at  St.  John's, 
Newfoundland  to  measure  the  "equivalent  slab  thickness"  of 
the  ionosphere  above  St.  John's.  Many  nights  during 
December,  1971  exhibited  appreciable  decreases  in  slab 
thickness  ('trsTEC/Nmax) . Specifically,  there  were  increases 
in  Nmax  on  these  nights  accompanied  by  constant  or 
decreas i ng  TEC. 

As  in  the  previous  chapter,  the  term  "TEC,"  will  refer 
to  the  equivalent  vertical  TEC  through  the  420-km  point, 
i.e.,  TEC  obtained  from  O and  the  constant  R value  at  420km 
height  along  the  ray  path.  Thus  the  ter. ns  "TEC"  and  "Faraday 
rotation"  (or  "Q"),  may  be  used  interchangeably. 

In  figure  IV-l  Nmax  and  TEC  have  been  plotted  vs. 
Universal  Time  for  five  nights  during  which  the  increasing 
Nmax-cons tant  TEC  effect  mentioned  above  occurred.  Note  that 
local  midnight  occurs  at  St.  John's  at  3:31  UT,  and  at  4:01 
UT  at  Goose  Bay.  The  Kp  Index  also  has  been  plotted.  Note 
that  with  the  exception  of  December  9,  the  Increase  in  Nmax 
is  not  associated  with  an  increase  in  magnetic  activity 
according  to  the  Kp  index.  The  AE  index  and  H at 
Narssarssuaq  (61*  2N,  43®  4W  ) are  plotted  in  addition  to  the 


-constant  TEC  effect  occurred 


1 > " 
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above  three  parameters  for  the  same  five  days  in  figures 

IV-2  to  IV-6.  H is  the  horizontal  component  of  the 

• 

geomagnetic  field.  The  downward  arrow  indicates  the 

direction  of  negative  deflection  of  H.  Some  additional  data 
are  also  plotted  In  the  graph  for  Dec.  9 and  Dec.  15  which 
we  will  turn  to  in  sec.  4.3.  We  can  notice  that  (again  with 
the  exception  of  Dec.  9)  the  Umax  increase  is  not  associated 


with  an 

increase  in 

the  AE  index. 

A substorm  is 

observed 

at 

Narssarssuaq  on  all 

five  nights,  but  except  for 

Dec.  9, 

the 

recorded 

substorm 

occurs  about 

three 

hours 

before 

the 

1 ncrease 

in  St.  John's  Umax.  We 

wi  1 1 

address 

this  topic 

fur  ther 

when  we  di 

scuss  the  mechanism 

for  producing 

the 

observed  slab  thickness  decreases  (sec.  4.4). 

Let  us  now  quantify  this  effect  for  December,  1971,  In 
section  4.2  the  results  of  a statistical  survey  of  increases 
in  foF2  are  presented.  For  December,  1971,  it  was  found  that 
at  least  a 0.  3MHz  increase  in  foF2  was  observed  on  21  of  25 
nights.  A comparable  sized  increase  in  TEC  is  about  20*. 
increases  of  more  than  20*  in  TEC  occurred  on  only  6 out  of 
27  nights  for  which  TEC  measurements  were  available.  To  get 
a feeling  for  the  variability  of  Umax  and  TEC  for  the  month 
as  a whole,  these  two  parameters  are  graphed  in  the  format 
of  an  overplot  In  figure  IV-7.  On  the  left  half  of  the 
figure,  the  two  parameters  are  plotted  for  all  of  the  days 
for  which  data  were  available.  To  the  right,  ten  of  the  days 
on  which  substantial  increases  In  Nmax  were  observed  are 
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Figure  I V- 5 Several  parameters  plotted  for  the  night  of 
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Overplots  of  Nmax  and  TEC  for  December,  1971 


plotted.  For  the  monthly  overplot,  the  data  collected  during 
the  great  magnetic  storm  of  Dec.  17  stand  out  sharply  from 
the  other  data  as  the  only  night  on  which  TEC  rose  above 
10l*el/cm*  and  Nmax  rose  above  3xl0Sel/cm*.  With  the 
exception  of  this  day,  an  "envelope"  drawn  around  the  plots 
would  extend  from  about  2.5  to  5.0  for  TEC  and  from  about  5 
to  15  for  Umax  (in  the  units  shown).  This  demonstrates  the 
significantly  greater  variability  of  peak  density  compared 
with  TEC  at  L=4  during  winter  nighttime  conditions. 
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4.2  Statistical  Survey  of  Winter 
Nighttime  Increases  in  foF2 

For  the  determination  of  slab  thickness,  measurements 
of  both  TEC  and  Nmax  are  required.  Since  the  TEC 
observations  from  Goose  Bay  commenced  in  November,  1971,  our 
knowledge  of  the  stab  thickness  parameter  at  St.  John's  is 
limited  to  the  period  after  that  month.  It  would  be  foolish, 
however,  to  ignore  the  wealth  of  foF2  data  available  before 
November,  1971,  in  our  attempt  to  understand  ionospheric 
processes  near  L»4.  We  have  gathered  24  years  of  St.  John's 
foF2  data  from  the  World  Data  Center  A for  Solar  Terrestrial 
Physics  and  from  the  Canadian  Dept,  of  Communi nations. 
Telecommunications  Engineering  Branch.  The  months  for  which 
data  were  analysed,  the  frequency  of  the  data  (whether  we 
obtained  foF2  values  at  15  minute,  half-hourly  or  hourly 
Intervals),  and  the  source  of  the  data  are  given  in  Table 
IV-  I . 

As  can  be  seen  from  Table  I V— I , the  December  data  we 
collected  spans  three  solar  cycles.  To  check  the  seasonal 
dependence  of  the  nocturnal  foF2  increases,  we  also  examined 
two  six  month  periods,  October,  19  71  to  March,  19  72,  and 
October,  1972  to  March,  1973.  The  results  of  this  survey 
provide  a quantitative  summary  of  the  magnitude  of  the  foF2 
increases,  as  well  as  their  seasonal  and  solar  cycle 
depe ndence. 
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Table  IV- 1 

Data  Used  In  Statistical  Survey 


Per i od 

Data  Interval 

Source 

of  Data 

Oct. 

19  71-March, 

1972 

15  minute 

C 

Oct. 

1972-March, 

1973 

15  minute 

C 

Dec. 

1945,  1947, 

and  1948 

hour  1 y 

w 

Dec. 

1950,  1951, 

and  1953 

hour  1 y 

w 

Dec. 

1957,  1958, 

and  1959 

hour  1 y 

w 

Dec. 

I960,  1901, 

and  1962 

hour  1 y 

w 

Dec. 

1963,  and  1964 

hour  1 y 

w 

Dec. 

1965  and  1966 

30  minute 

c 

Dec. 

1967,  1968, 

and  1969 

15  mi nute 

c 

Dec. 

1970,  1971, 

and  197/ 

15  minute 

c 

Dec. 

19  73 

15  mi nu  te 

c 

Dec. 

1974 

hour  1 y 

c 

Code 

for  sources 

of  data: 

C 3 Canadian  Dept,  of  Communl cations#  Telecommunications 
Engineering  Branch 


W * World  Data  Center  A for  Solar  Terrestrial  Physics 
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Each  night  was  examined  for  primary  and  secondary 
character  I stl cs.  Each  night  was  assigned  one  of  five  primary 
character  I stl cs  with  the  following  selection  criteria: 

low  bump:  Increase  from  .trough  to  crest  given  by: 
.3MHzi*foF2  1.5MHz 

bump:  Increase  from  trough  to  crest  given  by: 

. 6MHzi  Af  oF2  1.  9MHz 

large  bump:  increase  from  trough  to  crest  1 MHz  or 

greater . 

no  bump:  Increase  less  than  ,3MHz/  less  than  1 hour  long 
or  based  on  one  data  point/  except  for  cases  with  frequency 
of  data  half  hourly  or  hourly/  when  a single  data  point 
could  determine  the  primary  character  I stl c. 

Insufficient  data:  unable  to  determine  whether  an  Increase 
occurred. 

After  each  night  was  placed  In  one  of  the  above  five 
categories/  a secondary  characteristic  was  sometimes 
assigned.  In  determining  secondary  character  I s tl cs,  the 
criterion  that  It  must  be  possible  to  judiciously  eliminate 
any  single  data  point  without  eliminating  the  effect  was 
applied.  Failing  this  test*  the  secondary  characteristic  was 
not  assigned  except  when  only  hourly  or  half  hourly  values 
were  available.  The  secondary  characteristics  and  their 
selection  criteria  follow: 

Irregular  bump:  more  than  1 peak.  Trough  In  between  at 
least  0.3  MHz  below  the  two  peaks. 


short  bump:  “bump"  less  than  3 hours  long  (If  only  hourly 
values  available,  less  than  4 hours  long). 


plateau:  less  than  0.1  MHz  variation  for  at  least  3 hours 


at  the  maximum  of  the  Increase  (no  plateaus  when  only  hourly 


Pages  of  the  Canadian  data  are  reproduced  showing  each 


of  these  effects  in  figures  IV-8  to  IV-12.  Numerical  Tables 


of  the  seasonal  results  and  the  yearly  results  are  provided 
as  Table  I V— I i and  Table  I V— I V respectively. 

In  both  tables,  the  number  of  nights  with  each 
characteristic  is  listed  for  each  month.  The  last  column 
gives  the  fractional  occurrence  of  the  effect,  l.e.,  the 
number  of  nights  with  an  Increase  of  0.3  MHz  or  more 
divided  by  the  number  of  nights  with  sufficient  data  to 
determine  whether  an  foF2  increase  occurred.  For  Table 
IV-II,  the  numbers  for  October,  1971  and  October,  1972  have 
been  added  and  the  fractional  occurrence  determined  for  the 
composite  "mean  October."  The  same  has  been  done  for  all  the 
other  months  in  the  table.  The  mean  monthly  fractional 
occurrences  of  the  foF2  are  plotted  in  figure  I V- 13 . 

The  fractional  occurrence  in  December  Is  greatest, 
while  that  In  October  and  March  is  least.  The  sharp  dip  In 
January  is  probably  not  significant,  considering  the  fact 
that  only  two  years'  worth  of  data  were  used.  In  spite  of 
this  small  data  base,  we  can  initially  conclude  that  the 
effect  we  have  been  studying  is  a winter  one,  with  the 
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Table  I V- 1 1 

Statistical  Survey  of  Nocturnal  Winter  foF2  Increases 

SEASONAL  ANALYSIS 
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greatest  fractional  occurrence  most  likely  to  be  found  in 
December. 

A convenient  tool  for  quantifying  the  correlation 
between  two  or  more  variables  is  the  Pearson  Product-Moment 
correlation  coefficient.  The  correlation  between  a variable 
and  the  same  variable  using  this  method  is  1#  while  the 
correlation  between  a variable  and  its  inverse  (x  and  1/x) 
would  be  -1.  If  two  variables  are  independent  of  one  another 
such  that  the.  e is  no  correlation#  the  correlation 
coefficient  Is  zerc- 

This  method  has  been  used  to  determine  the  correlation 
between  solar  zenith  angle  (angular  distance  of  the  Sun 
from  the  zenith)  at  noon  and  the  primary  and  secondary 
characteristics  (with  the  exception  of  plateaus#  because  of 
the  small  number  of  them).  The  results  are  given  below  in 
Table  I V- I 1 1 , 

Table  I V- I I I 

The  correlation  between  noontime  solar  zenith  angle  and: 
fractional  occurrence  of  all  Increases  * .#8 
fractional  occurrence  of  low  bumps  * .46 
fractional  occurrence  of  bumps  « .91 
fractional  occurrence  of  large  bumps  * .83 
fractional  occurrence  of  Irregular  bumps  e .27 
fractional  occurrence  of  short  bumps  = -.21 


The  high  correlation  for  fractional  occurrence  of  all 
Increases  confirms  our  Initial  conclusion  that  the  effect  we 
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have  been  studying  is  a winter  one.  There  seems  to  be  a 
significant  difference  between  the  correlation  coefficients 
for  low  bumps  (.46)  and  for  the  larger  bumps  (.91  and  .83). 
One  can  notice  from  Table  IV-II  that  a large  number  of  low 
bumps  occurred  in  October.  February,  and  March,  while 
relatively  few  large  bumps  occurred  during  those  months.  The 
correlation  coefficients  for  the  secondary  character! s ti cs 
are  small,  and  since  they  represent  a smaller  fraction  of 
the  nights,  we  must  discount  their  importance.  We  will 
suggest  a possible  reason  for  this  difference  between  the 
correlation  coefficients  for  the  low  and  larger  bumps  in 
section  4.4.  when  we  turn  our  attention  toward  an 
explanation  of  the  mechanism  for  the  foF2  increases. 

Twenty  four  years  of  December  statistics  are  compiled 
in  Table  IV-IV.  In  addition  to  appearing  in  Table  IV-IV.  the 
fractional  occurrence  of  all  increases  0.3  MHz  and  greater 
is  plotted  in  figure  IV-14.  A glance  at  this  figure  does  not 
seem  to  reveal  any  correlation  with  the  solar  cycle. 

As  was  done  for  the  seasonal  case,  correlation 
coefficients  were  computed,  this  time  between  observed 
sunspot  number  and  the  primary  and  secondary  characteristics 
(again  excluding  plateaus).  The  results  are  given  >n  Table 
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Statistical  Survey  of  Nocturnal  Winter  foF2  Increases 


SOLAR  CYCLE  DEPENDENCE 


Dec.  1945  749 


Dec.  1957  4 5 15  5 
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Table  IV-IV  (continued) 
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Table  IV-V 

The  correlation  between  observed  sunspot  number  and: 

fractional  occurrence  of  all  increases  ■ -.13 

' 

fractional  occurrence  of  low  bumps  « -.49 
fractional  occurrence  of  bumps  * -.23 
fractional  occurrence  of  large  bumps  * .48 
fractional  occurrence  of  irregular  bumps  « -.12 
fractional  occurrence  of  short  bumps  » .26 

There  is  no  significant  correlation  between  sunspot 
number  and  fractional  occurrence  of  all  increases.  One 
striking  difference  is  the  correlation  coefficients  for  low 
bumps  and  for  large  bumps,  the  largest  two  of  all  the  above 
coefficients  and  strikingly  opposite  one  another.  When  these 
two  variables,  fractional  occurrence  of  low  bumps  and 
fractional  occurrence  of  large  bumps  are  correlated  with 
each  other,  the  result  Is  a correlation  coefficient  of  -.70, 
a significant  result. 

As  in  our  sunspot  cycle  correlation  study,  the 
correlation  coefficients  for  the  secondary  characteristics 


do  not 

appear 

to  be  significant. 

based 

as 

they 

are 

on  a 

smal ler 

amount 

of  data. 

We 

have. 

in  conclusion,  examined 

the 

seasonal 

and 

sunspot 

cycle 

dependence  of  the 

foF2 

I ncr eases 

with 

the 

following  results:  there  Is  a strong  preference 

for 

this 

effect  to  occur  In  winter  (November-February ) but  no  way  of 
predicting  whether  a given  year  will  have  a large  fractional 


4 


101 


< 

I 


occurrence  of  the  overall  effect;  there  are  significant 
differences  between  the  seasonal  and  sunspot  cycle 
dependences  of  low  and  large  bumps,  the  large  bumps  more 
strongly  clustered  in  mid-winter  and  more  likely  to  occur 
during  high  sunspot  years. 
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4.3  Case  Studi es: 

December  9 and  December  15,  19  71 

4,3.1  Description  of  Available  Data 

Two  nights,  December  8-9,  1971,  and  December  14-15, 
1971,  were  chosen  for  more  detailed  analysis  because  of  the 
relatively  large  variety  and  quantity  of  geophysical  data 
available  for  these  two  nights,  in  addition  to  information 
plotted  in  figures  IV-3  (for  Dec,  9)  and  I V— 6 (for  Dec.  15), 
Isis  II  topside  sounder  passes  recorded  Ne(h)  vs.  latitude 
on  both  nights.  Two  passes  are  of  value  for  our  purposes  on 
December  9.  The  first  occurred  at  about  4:20  UT  and  at 
longitude  71  degrees  west,  while  the  second  occurred  at 
about  6:15  UT  and  about  100  degrees  west.  Two  passes  are 
also  of  interest  on  December  15,  occurring  at  about  2:25  UT 
and  49  degrees  west  longitude,  and  at  about  4:20  UT  and  73 
degrees  west  longitude. 

While  discussing  the  data  available  for  our  case 
studies,  we  should  take  a moment  to  consider  the  question  of 
the  quality  of  these  data.  As  mentioned  in  section  2. 1.3. 3, 
the  quality  of  the  Isis  II  lonograms  varied  greatly,  with 
real  height  obtained  from  the  best  quality  lonograms 
accurate  to  about  *5 %. 

The  primary  evidence  for  the  foF2  increases  Is  provided 
by  lonograms  taken  at  St,  John's,  Newfoundland.  We  obtained 
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foF2  values  for  many  of  the  same  ionograms  from  different 
sources  In  order  to  check  the  accuracy  In  foF2.  The  three 
sources  compared  were:  the  Canadian  source  previously  cl  ted. 
World  Data  Center  A,  previously  cited,  and  Dr.  J.  Buchau  at 
the  Air  Force  Cambridge  Research  Laboratories.  The  foF2 
scaled  by  the  three  sources  always  agreed  to  within  0.2  MHz. 

Two  ionograms  taken  during  the  Dec.  13  event  Cone 
before,  and  one  after  the  greatest  part  of  the  increase 
occurred)  are  reproduced  In  figure  IV-15.  While  the  two 
ionograms  are  not  of  the  best  quality,  a sizeable  Increase 
in  foF2  Is  apparent  from  the  first  ionogram  to  the  second. 
The  foF2  values  marked  on  the  figure  are  those  scaled  by  Dr. 
J.  Buchau. 

Let  us  now  turn  our  attention  to  figure  I V—  3 for 
December  9,  1971.  As  previously  noted,  figure  I V- 3 shows 
that  on  Dec.  9,  the  most  disturbed  magnetic  conditions 
coincided  with  the  onset  of  the  increase  In  peak  density.  A 
much  larger  increase  began  about  an  hour  later  at  Maynard, 
which  is  18.8  degrees  west  of  St.  John's.  The  Increase  in 
TEC  seen  at  Sagamore  Hill  should  be  compared  with  Nmax  at 
Wallop's  Island,  near  the  Sagamore  Hill  420-km  point,  rather 
than  with  Nmax  at  Maynard.  Unfortunately,  no  data  were 
available  for  the  time  in  question  at  Wallop's  Island. 

On  the  St.  John's  Nmax  vs.  UT  plot,  two  times  are 
marked  (3:30)  and  P*  (6:00).  A least  squares  fit  between 
the  profile  calculated  by  World  Data  Center  A and  the 
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semi-parabola  used  In  the  computer  model  gave  a best 
yzero=213km  at  P|  and  164km  at  P4.  One  can  see  from  the 
plots  of  these  two  profiles  in  figure  IV-3,  that  the 
bottomside  has  filled  In  considerably/  which  seems  to  Imply 
a downward  flux  of  plasma. 

On  the  night  of  December  8-9/  no  Ne(h)  vs.  latitude 
Isis  II  data  were  collected  near  the  longitude  of  St. 
John's.  We  do,  however,  have  the  /I  degree  west  longitude 
pass  previously  mentioned,  Nmax,  N(700),  and  N(10U0)  vs. 
geographic  latitude  from  this  pass  are  plotted  In  figure 
IV-16.  Slight  changes  were  made  at  two  points  on  the  array 
where  Nmax  and  hmax  values  were  derived  from  ground  based 
lonosonde  data.  Note  the  extremely  sharp  poleward  trough 
edge.  There  is  an  unfortunate  gap  in  the  data  between  51.68 
degrees  and  66.56  degrees.  The  computer  model  Interpolation 
between  these  two  points  might  exaggerate  the  actual  plasma 
content  In  this  latitude  range. 

Let  us  now  discuss  the  December  15  data. 

In  figure  I V- 6 we  can  note  that  the  beginning  of  the 
Nmax  Increase  coincided  (as  on  the  case  of  December  9)  with 
the  greatest  magnetic  activity,  A much  larger  Nmax  Increase 
did  not  begin,  however,  until  about  three  hours  later.  The 
Nmax  vs.  UT  plot  for  Wallop's  Island  looks  similar  to  that 
for  St.  John's.  Unlike  the  case  of  Dec,  9,  when  the  increase 
at  Maynard  was  far  greater  than  that  at  St,  Johns,  on  the 
15th,  the  Increase  to  the  south  (at  Wallop's  Island)  is  much 
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less  than  that  at  St*  John's.  The  increase  at  Wallop's 
Island  occurs  later  than  the  Increase  at  St.  John's,  but  at 
about  the  same  local  time. 

Profiles  for  the  times  3:00  UT,  4:15  UT,  and  6:45  UT 
obtained  from  World  Data  Center  A are  shown  and  Indicated  on 
the  Nmax  vs.  UT  plot  for  St.  John's  as  P|,  Pa,  and  Pj  . A 
least  squares  fit  between  the  profile  calculated  from  the 
St.  Johns's  lonogram  by  World  Data  Center  A and  the  best  fit 
semi-parabola  used  in  the  computer  model  gave  yzero*219km  at 
Pt,  yzero»153km  at  P*,  and  yzero»2U3km  at  Pj  . Significant 
changes  occurred  on  the  bottoms ide  during  this  event  as 
these  results  Indicate. 

It  is  very  fortunate  that  we  have  available  Isis  II 
data  near  the  longitude  of  St.  John's  on  December  15.  Data 
from  this  2:25  UT  pass  at  49  degrees  west  geographic 
longitude  are  plotted  in  figure  I V- 17.  One  change  in  this 
array  was  made,  viz.,  Nmax  and  hmax  obtained  from  St.  John's 
ionograms  were  substituted  for  the  topside  values. 

4.3.2  Computer  Simulations  of  Slab  Thickness 
Decreases  for  Dec.  9 and  Dec.  15,  1971 

4, 3. 2.1  Preliminary  Considerations 

As  has  been  seen,  the  computer  model  we  have  developed 
takes  as  Input  data  an  Ne(h)  vs.  latitude  array,  and 
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computes  the  Faraday  rotation  from  the  geostationary 
satellite  to  a given  station.  We  wish  to  use  this  model  to 
reproduce  the  effect  observed  on  Dec.  9,  and  Dec.  15,  i.e., 
constant  Faraday  Rotation  accompanied  by  increasing  peak 
density  at  St.  John's. 

Although  we  do  have  an  Ne(h)  vs.  latitude  array  near 
the  longitude  of  St.  John's  for  Dec.  15,  the  array  for  Dec. 
9 closest  in  longitude  to  St.  John's  is  at  71.3  degrees 
west,  too  far  to  the  west  for  us  to  expect  to  reproduce  the 
St.  John's  and  Goose  bay  results  very  accurately.  It  is 
probably  reasonable  to  assume,  however,  that  a similar  array 
would  be  found  from  a pass  near  St.  John's.  Therefore,  it 
would  be  of  gr?at  value  for  us  to  try  to  reproduce  the 
observed  slab  thickness  decrease  which  occurred  on  Dec.  9 as 
well  as  on  Dec.  15. 

We  might  ask,  at  this  point,  how  well  the  total 
electron  content  of  a vertical  column  at  420km  height  along 
the  ray  path  calculated  by  our  computer  model,  agrees  with 
the  actual  equivalent  vertical  TEC  observed  at  Goose  Bay.  In 
Table  IV-VI  is  listed  observed  and  calculated  equivalent 
vertical  TEC  from  Goose  Bay  and  Hamilton,  Mass.  (Sagamore 
Hill  Radio  Observatory).  Notice  that  the  discrepancy  between 
observed  and  calculated  TEC  is  larger  at  Goose  Bay  than  at 
Hamilton  for  both  dates.  This  may  be  due  to  geometry  effects 
not  fully  taken  into  consideration  in  our  model.  The 
computer  simulation  assumes  that  the  satellite  is  on  the 


Table  IV-VI 

Calculated  and  Observed  TEC 

Equivalent  Vertical  TEC 
xlO'^el/ cm 1 


Date 

Statl on 

Observed 

Time  of 

Calc. 

Time 

of 

TEC 

Obser vatl on 

TEC 

Isis 

1 1 

pass 

Dec. 

9 Goose  Bay 

3.5 

4:00  UT 

2.43 

4:20 

UT 

(53°.  3N,  60°.  3W) 

3,5 

5:00  UT 

Hami  1 ton 

4.7 

4:00  UT 

3.95 

(42°.  6N,  70°.  8W  ) 

4.7 

5:00  UT 

Dec. 

15  Goose  Bay 

4.4 

2:00  UT 

3.04 

2:25 

UT 

(53°.  3N,  60°.  3W) 

4.3 

3:00  UT 

Hami 1 ton 

6.3 

2:00  UT 

5.54 

(4  2°.  6N,  70°,  8W) 

5.7 

3:00  UT 

Ill 

station's  meridian.  For  ATS-3,  stationed  at  7d°W,  this  is  a 
good  approx! mat! on  for  the  Sagamore  Hill  case/  but  not  for 
the  Goose  Bay  ray  path. 

An  additional  factor  contributing  to  the  remaining 
discrepancy  between  observed  and  calculated  equivalent 
vertical  TEC  is  found  by  remembering  that  the  lowest  height 
at  which  electron  density  is  given  by  the  Isis  II  topside 
sounder  is  often  somewhat  above  hmax  (see  sec.  2. 1.3. 3). 
However,  the  density  at  this  height,  which  is  lower  than 
Nmax,  is  taken  in  our  model  to  be  the  peak  density.  This  may 
represent  a source  of  systematic  error,  suotractlng  a 
considerable  amount  of  plasma  around  the  F-2  peak  in  our 
model  and  lowering  the  calculated  Faraday  rotation. 

Of  course,  we  cannot  expect  perfect  agreement  oetween 
observed  and  calculated  TEC  for  other  reasons.  The  computer 
model  reads  in  only  five  parameters  to  represent  the  entire 
ionosphere.  The  exospheric  models  are  not  perfect,  either, 
extrapolating  as  they  do  from  the  very  idealized  conditions 
at  lOUUkm  all  the  way  to  the  geostationary  satellite.  In. 
addition,  whenever  the  satellite  pass  is  at  a different 
longitude  from  the  ray  path  to  the  geostationary  satellite, 
different  conditions  are  being  sampled. 

In  spite  of  these  sources  of  error,  we  can  be  confident 
that  the  semi -empi r i ca I computer  model  we  have  developed 
should  enable  us  to  reproduce,  at  least  qual  I tati  vely, 
observed  physical  effects  which  occur  on  a large  scale. 
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We  will  use  this  computer  model  in  two  ways  to  attempt 
to  reproduce  for  Dec.  9 and  Dec.  15,  the  observed 

anomalistic  effect,  i.e.,  constant  TEC  observed  at  Goose  Bay 
coincident  with  an  increase  in  peak  density  at  St.  John's. 
The  changes  which  must  be  accounted  for  are  shown  in  Table 
IV-VI I and  Table  I V-V I I I. 

The  first  method  which  we  will  use  involves  a 
latitudinal  shift  of  the  entire  Ne(h)  vs.  latitude  array  (as 
was  done  in  Chapter  2).  Shifting  the  poleward  edge  of  the 
electron  density  trough  over  St.  John's  will  clearly 
Increase  Nmax,  We  must  discover  the  effect  of  such  a shift 
on  the  Faraday  rotation. 

The  second  method  Is  a vertical  distortion  of  the  Ne(h) 
vs.  latitude  array  without  a latitude  shift  of  the  array  as 
a whole.  It  appears  fairly  likely  that  there  are  a large 
number  of  latitude  dependent  distortion  schemes  which  will 

reproduce  the  desired  effect.  One  such  scheme  will  be 

* 

presented  which  seems  physically  reasonable. 

After  using  both  methods,  results  from  a certain 
combination  of  the  two  \/i  1 1 also  be  presented. 


4. 3. 2. 2 Latitude  Shift  Method 


The  reader  will  recall  that  In  Chapter  2,  in  order  to 
examine  the  effect  of  a severe  trough  on  the  determi nation 
of  an  appropriate  M factor,  the  feature  of  our  computer 
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4:20  UT 

6:00  UT 

% change 

Umax  5.4  C 1 s I s II) 

7.1 

♦ 31 

xl0#el/cml 

TEC  3.5 

3.4 

-3 

xlO14  el/cm* 

Table  I V-V I 1 1 

OBSERVED  INCREASES  - December  15 


2:25  UT 

6:45  UT 

4 change 

Umax  9.1  (Isis  II) 

17.9 

♦97 

xl0*el/cm* 

TEC  4.4 

3.65 

-17 

xl0|a,el  /cm* 

114 


model  which  allows  us  to  shift  the  Ne(h)  vs.  latitude  array 
an  arbitrary  number  of  degrees  north  or  south  was  employed. 
With  this  technique  available  to  us,  we  now  attempt  to 
reproduce  the  observed  temporal  behavior  of  Nmax  and  Cl  on 
Dec.  9 and  Dec.  Is,  1J71  by  shifting  latitudinal  Ne(h) 
arrays  southward.  The  results  of  shifting  the  Dec  9th  array 
5 times  southward  In  one  degree  steps  are  shown  in  Table 
IV-IX.  For  yzero,  we  chose  the  value  of  21akm  for  all 
latitudes  In  the  input  array  except  the  two  northernmost 
points,  which  are  north  of  the  trough  poleward  edge.  There  a 
value  of  165km  was  used.  The  motivation  for  these  choices  of 


yzero 

is  apparent  from 

remember i ng 

that  yzero 

for  Pt 

i n 

f l gure 

IV- 3 was  213km 

and  yzero 

for  Pa  was 

164km. 

If 

shifting  the  array  south  is  to  represent  a temporal  process, 
then  yzero  at  St.  John's  must  change  from  about  215km  to 
about  165km  as  the  poleward  edge  passes  over  St.  John's. 

One  can  see  that  In  Table  IV-IX,  changing  from  a shift 
of  -3  degrees  to  one  of  -5  degrees  (which  corresponds  to 
moving  the  poleward  edge  of  the  trough  over  St.  John's), 
gives  a 395£  increase  in  Umax  coincident  with  a 29*  increase 
In  While  this  Faraday  rotation  increase  is  very  small 
compared  to  the  Umax  Increase,  It  is  still  significant,  and 
inconsistent  with  observations.  Fu rthermore,  to  reproduce 
the  observations  we  must  compare  a zero  degree  shift  with 
another  shift,  assumed  to  give  the  ionospheric  densities  at 
a later  time.  Examination  of  Table  IV-IX  shows  that  the 
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Table  IV- IX 

LATITUDE  SHIFT  RESULTS  - December  9 


yzero“165km 

north  of  poleward 

edge. 

yzero«215km 

el  sewhere. 

Shi  ft 

St.  John's  Nmax 

Goose  Day  O 

TEC 

xlO^el/cm1 

degrees 

xlO1*  el /cm* 

0 

5.5 

100 

2.4 

-1 

5.3 

127 

2.9 

-2 

5.  3 

204 

4.6 

-3 

6.3 

277 

6.2 

-4 

13.6 

325 

' 7.3 

-5 

31.2 

357 

3.0 

Tab  1 e 

i IV-X 

LATITUDE  SHIFT  RESULTS 

- December  9 

yzero=230km 

north  of  poleward 

edge. 

yzero*120  km 

In  trough 

south  to  St.  John's,  yzeroa215km  from  St.  John' 

's  southward 

Shift 

St.  John's  Nmax 
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xl04e l /cm* 

degrees 

x 1 0 IZ  e 1 / cm' 

0 

5.5 

130 

2.9 

-1 

5.3 

ItO 

3.2 

-2 

5,3 
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3.7 

-3 

6.  3 

219 

4.9 

-4 

19.2 

263 

5.9 

-5 

31,2 

291 

6.5 
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desired  effect  has  not  been  reproduced,  since  initially  (as 
we  shift  one,  two,  three  degrees  south),  Nmax  does  not 
increase  while  Cl  does,  contrary  to  the  actual  temporal 
sequence  of  events.  In  addition,  comparing  the  zero  degree 
shift  with  the  -5  degree  shift  gives  a 467%  Nmax  increase 
and  a 228%  Cl  increase,  completely  at  variance  with  the 
observed  increases  listed  in  Table  I V- VII. 

If  we  suppose  that  the  change  in  yzero  at  the  trough 
poleward  edge  is  very  sharp,  taking  yzero*230km  north  of  it 
and  12dkm  south  of  I t as  far  as  St.  John's  (with  yzero*215km 
southward  from  St.  John's),  then  we  can  reduce  the  Faraday 
rotation  increase  in  changing  from  no  shift  to  a -5  degree 
shift  by  about  50%.  This  Is  because  although  Nmax  in  the 
trough  is  very  low,  yzero  is  required  to  be  very  low  there 
so  that  the  vertical  TEC  is  greater  than  it  otherwise  would 
be.  Even  with  this  artificial  device  we  cannot  reproduce  the 
observations  for  we  still  have  a 12**%  increase  in  Faraday 
rotation  coupled  with  the  467%  Nmax  increase.  The  results 
for  the  shifts  from  0 to  -5  degrees  using  these  values  for 
yzero  are  listed  in  Table  IV-X. 

For  the  Dec.  15  case  we  initially  tried  yzero*2U3km 
north  of  the  poleward  edge,  153km  south  of  the  poleward  edge 
to  the  latitude  of  St.  John's  and  222km  from  St.  John's 
southward.  As  we  shift  southward,  yzero  above  St.  John's 
changes  from  2 2 2 kin  to  153km  to  203km,  which  Is  the  observed 
temporal  behavior. 
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The  results  for  shifts  of  zero  to  -6  degrees  In  one 
degree  steps  for  Dec.  15  are  shown  In  Table  I V— X I • A shift 
of  six  degrees  Is  necessary  for  the  poleward  edge  to  reach 
St.  John's  producing  a 724  Increase  In  Nmax.  Faraday 
rotation  increased  by  a much  smaller  percentage/  234. 

We  now  artificially  distort  the  trough  bottomslde 
profiles  north  of  St,  John's  by  letting  yzero  equal  120km 
there  Instead  of  153km.  This  decreases  the  Faraday  rotation 
Increase  from  the  zero  degree  shift  to  the  -6  degree  shift 
to  174  (see  Table  I V—  XII).  A further  distortion  Is  obtained 
by  keeping  yzero«120km  south  of  the  poleward  edge  to  St. 
John's  and  letting  yzerom230km  north  of  the  poleward  edge 
instead  of  203km,  Now  the  724  Nmax  Increase  from  0 degree 
shift  to  -6  degree  shift  Is  accompanied  by  a 7t  Faraday 
rotation  Increase  (see  Table  IV-XIII).  A summary  of  these 
three  sets  of  results  Is  found  In  Table  I V—  XIV.  Graphs  of 
Faraday  rotation  and  Nmax  for  the  three  calculations  are 
shown  In  figure  IV-18  for  a wide  range  of  latitude  shifts. 
The  best  agreement  to  the  observeJ  effect  on  Dec.  15  (small 
TEC  decrease  coincident  with  Nmax  Increase)  Is  found  In  the 
third  calculation  (74  CX  Increase). 

To  get  a feeling  for  the  change  In  the  latitudinal 
structure  sampled  by  the  ray  path/  two  contour  plots  are 
drawn  in  figure  IV- 19,  one  for  the  zero  degree  shift  and  one 
for  the  -6  degree  shift,  for  the  results  of  the 
aforementioned  third  calculation.  The  ray  path  is  drawn 
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Table  IV-XI 

LATITUDE  SHIFT  RESULTS  -December  15 

yzero*203km  north  of  poleward  edge,  yzero»153km  In  trough 
south  to  St.  John's,  yzero*222km  from  St.  John's  southward. 
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LATITUDE  SHIFT  RESULTS  - December 
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yzero* 

Table  IV-XI 1 1 

Latitude  Shift  Results  - December  15 
230km  north  of  poleward  edge,  yzero*120km 
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Figure  IV-18  Model  Nmax  and  Cl  for  Dec,  15  using 
various  yzeros  for  a wide  range  of 
latitude  shifts 


yzero*203km  north#  153km  south  of  poleward  edge 

yzero*203km  north#  120km  south  of  poleward  edge 

yzero“230km  north#  120km  south  of  poleward  edge 
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through  the  interpolated  density  profiles  of  Ne(h)  vs. 
latitude  from  0 km  In  height  to  lJOOkm  in  height,  which  were 
provided  by  the  computer  model.  Notice  that  the  poleward 
edge,  at  the  extreme  rightCnorth)  in  the  zero  degree  shift 
plot,  has  been  moved  to  cross  the  ray  path  near  400km  height 
which  greatly  increases  the  Faraday  rotation  contribution 
there.  This  is  compensated  for  to  a great  extent  by  two 
factors:  the  decrease  in  electron  densities  at  greater 
heights  along  the  ray  path  due  to  the  fact  that  this  is  now 
the  trough  region;  and  the  higher  yzeros  in  the  poleward 
region  where  the  ray  path  is  traveling  through  lower 
hei ghts . 

A glance  at  Table  IV-VIII  will  reveal  that  while  we 
have  greatly  reduced  the  Faraday  rotation  increase,  we  still 
have  not  simulated  the  effect.  A 11%  decrease  in  Faraday 
rotation  cannot  be  duplicated  by  the  latitude  shift  method 
using  physically  reasonable  values  of  yzero.  The  values  of 
yzero  we  had  to  choose  to  cut  down  the  Faraday  rotation  were 
quite  arbitrary.  The  value  of  IzOkm  was  quite  a bit  lower 
than  the  values  which  we  have  said  were  observed  at  St. 
John's,  and  230km  is  higher  than  yzero  observed  at  6:45  UT. 

We  can  conclude  our  discussion  of  this  latitude  shift 
method  by  saying  that  equatorward  motions  of  latitudinal 
Ne(h)  profiles  cannot,  by  themselves,  explain  the  observed 
anomalistic  effect,  l.e,,  flat  or  decreasing  TEC  combined 
with  an  Nmax  increase.  This  conclusion  is  strengthened  by 
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examination  of  the  succeeding  Isis  II  passes  on  both  Dec.  9 
and  Dec.  15  (compare  figures  IV-20  and  IV-21  with  figures 
IV-16  and  I V— 17).  Although  the  succeeding  passes  are  far  to 
the  west  of  the  original  ones,  if  any  substantial  motion  of 
the  Ionization  trough  occurred.  It  should  show  up  in  the 
later  passes.  Examination  shows  that  the  poleward  edge  did 
not  seem  to  change  its  position  appreciably  on  Dec.  9,  and 
on  Dec.  15,  it  is  impossible  to  tell  where  the  poleward  edge 
was  at  the  later  time  (if  it  actually  was  present). 


In  spite 

of  this 

negat i ve 

conclusion 

i,  we  note 

that 

appreci able 

deer  eases 

in  slab 

thi ckness 

(TEC/Nmax) 

were 

reproduced 

on 

both  Dec 

. 9 and 

Dec.  15. 

Even  though 

TEC 

increased,  the  increase  in  Nmax  was  much  greater. 

In  short,  it  is  certainly  possible  that  some 
latitudinal  plasma  motions  occurred,  as  a simple  drift  model 
to  be  discussed  presently  suggests.  These  latitudinal 
motions  are  apparently  not  the  primary  cause  of  the 
phenomenon  we  are  concerned  with  (flat  TEC  with  increasing 
Nmax),  however. 


4. 3. 2. 3 Vertical  Drift  Method 

It  seems  reasonable  to  hypothesize  that  a vertical 
redistribution  of  plasma  over  a wide  latitude  range, 
increasing  Nmax  but  depleting  the  electron  density  by 
appropriate  amounts  at  certain  other  heights,  will  produce  a 
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constant  fi  (and  Its  derived  TEC)  coincident  with  the  N max 
Increase,  There  are  literally  an  Infinite  number  of  ways  to 
do  this.  The  Nmax  Increase  need  only  occur  over  a small 
latitude  range  to  reproduce  our  single  ground  based 
lonosonde  observations.  It  need  then  be  accompanied  by  a 
small  decrease  elsewhere  I n Ne  at  lower  and  higher  heights. 
If  one  increases  Nmax  over  a wide  latitude  range#  one  is 
still  at  liberty  to  specify  how  the  decrease  in  plasma 
density  at  other  heights  needed  to  keep  the  TEC  constant  Is 
to  be  distributed  along  the  ray  path. 

Even  though  succeeding  Isis  II  passes  are  about  30 
degrees  of  longitude  apart#  we  can  get  an  idea  of  whether 
vertical  distortion  occurred  over  a wide  latitude  range  by 
examining  N(700)  and  N(1000)  for  successive  passes.  One  can 
notice  by  comparing  figure  I V—  2 0 with  figure  IV-16  for  Dec. 
9#  and  figure  IV-21  with  figure  IV-17  for  Dec.  IS#  that  both 
N( 700 ) and  N(1000)  have  decreased  significantly  over  nearly 
the  entire  latitude  range. 

This  tells  us#  that  in  simulating  the  constant 
TEC- I ncr eas 1 ng  Nmax  effect  using  our  computer  model,  we 
should  distort  the  input  array  over  the  entire  latitude 


range. 
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A simple  model  for  "Magnetos pher i c Convection  at 
Ionospheric  Heights"  has  been  offered  by  Mend! llo  (1973).  In 
that  analysis  the  vertical  drift  Is  expressed  as  follows: 


In  this  formula: 

Et»  Is  a uniform  dawn  to  dusk  magnetos pher 1 c electric  field 
in  the  equatorial  plane. 

LT  is  the  local  time. 

tie£  is  the  magnetic  field  at  the  Earth's  surface  at  the  equator, 
h is  the  height. 
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change  in  NeCh)  are  related.  While  the  actual  relation  Is 
probably  quite  complicated/  the  following  simple  argument 
will  be  sufficient  to  provide  some  physical  justification 
for  the  particular  distortion  method  we  chose  for  our 
computer  simulation.  The  purpose  of  this  vertical  distortion 
exercize  Is  to  qualitatively  discuss  the  large  scale  changes 
of  ionospheric  structure  that  must  be  primarily  responsible 
for  the  observed  effects.  It  is  unfortunately  beyond  the 
scope  of  this  thesis  to  attack  the  full  time  and  spatial 
dependence  of  the  modified  ionospheric  structure. 

Consider  an  Imaginary  box  in  the  ionosphere  with 
cross-sectional  area  A and  height  Ah.  For  a downward  drift, 
the  fluxes  In  and  out  of  the  box  are  given  by: 


F; 


in 


(4.3) 


Font  * Vb.t 


Here  we  assume  no  horizontal  drifting  of  plasma,  and  for 
further  simplicity,  no  net  production  or  recombination  of 
free  electrons  In  the  box.  Then  we  can  write: 

Nk„  ( t+At ) AAh  * ( t ) AAh  ♦ ( Nt,f  VTtf  -Nk#t  ) AAt  (4.4) 

so  that: 

Nkw  ( t+At ) = Nu,  ( t ) ♦ ( NUf  VCop  Vk.t  )At/Ah  (4.5) 

Since  this  is  not  a large  box,  we  can  take 

vtof  “ vfc«C  a Wvert  (4.6) 


Then : 


^(t+At)  = Nu<t)  ♦ (AN^/Ah ) WvertAt  (4,7) 
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Here  AN^/Ah  should  be  an  average  value,  l.e.,  the  true 
value  of  AN/Ah  at  some  time  between  t and  t*&t. 
Unfortunately,  all  we  have  available  to  us  are  the 
Interpolated  density  arrays  for  the  beginning  of  the  effect 
we  wish  to  simulate.  Therefore  we  must  use: 

N(t+At)  * M(t)  ♦ (AN/Ah  )^Wver tAt  (4.8) 

A more  general  expression  for  dN/dt  is  given  by: 

-^  = c| - L + v- v N + Nv-v  (4.°)) 

This  is  the  continuity  equation  (1.10)  with  a motion  term 
added.  Here  we  took  q = L,  and  V*V=0,  so  (4.9)  reduces  to 
(4,8).  This  equation  should  be  applied  to  M(lOOO),  N(700) 
and  N max. 

It  is  clear  why  decreases  occur  at  7o0km  and  1000km, 
i.e.,  dN/dh  < 0.  To  see  why  Umax  increases,  where  dN/dh  = 0, 
consider  the  following,  We  cannot  apply  (4.8)  to  Nmax.  This 
is  because  dN/dh  = 0 at  hmax  for  any  profile,  and  in 

particular,  the  semi -par abol a we  have  modeled  around  the  F-2 
peak  is  such  that: 

N(hmax+Ah)  * N(hmax-Ah)  (4.1u) 

So  that  AN/Ah  is  zero  at  the  peak.  If  our  seni -parabol i c 
model  represented  the  correct  distribution,  we  could  not, 
according  to  (4.8),  Increase  Nmax  using  vertical  plasma 
drifts  (because  as  much  plasma  would  be  entering  the  region 
of  peak  density  as  would  he  leaving  it).  The  semi-parabola 
is  only  an  approximation,  of  course.  We  know  from  experience 
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that  generally  during  the  nighttime: 

N(hmax-Ah)  < N(hmax*Ah) 


Figure 

Nma* 

Nihmax  -aVi) 


We  have  drawn  a typical  nighttime  Ionospheric  profile 
schematically  in  figure  I V- 2 2 which  demonstrates  that  Ne(h) 
normally  falls  off  more  rapidly  below  hmax  (due  to  Increased 
recombination  rates)  than  it  does  above  hmax  (also  see 
figures  I V- 2 5 and  IV-26).  Therefore,  in  realityAN/Ah  Is 
positive  near  the  F-2  peak  and  according  to  (4.8),  Nmax  will 
increase  as  expected.  This  can  also  be  easily  seen  from 
(4.4)  with  < Nt<^.  Conversely,  A^/Ah  is  negative  at 

h»700km  and  h«l000km,  so  the  electron  density  decreases 
there,  as  desired. 

For  the  topside  heights,  we  thus  have: 

>'1iom(  t+At ) * N,<*.(t)  ♦ (AN/AhW  Wver  t(  1000  )At 
Nw  ( t*A*  > • N7oo<t>  ♦ (aN/a^W  Wvert(70Q)At  (4.11) 
Since  (4.2)  does  not  give  the  absolute  magnitude  of  Wvert, 
we  replace  WvertAt  In  (4.11)  by: 

K*w  * WvertAt  (4.12) 


where  from  (4.1): 

K=  ~EoDCos(SLT)At 

Be's  (4.13) 

K(E0Di  LT,  Be,s,&t) 

Then  equations  (4.11)  become: 

HtMoLt+At)  • Nm(t)  * (AN/Ahj>„,jtK*w(1000) 

Nioo(t^At)  * N7ao(t)  ♦ (aN/aI»^  K*w(700)  (4.14) 

The  equivalent  formalism  for  Nmax  would  be  an  equation 
of  the  form: 

Nmax  (t+ at)  = Nmax(t)  + (AN/ M^^h*^ax,tK*  *^(hmaO  (4|5) 

Nmaxit  + *t)  = Nmaxlt)  ^ | + w(UaO) 


Nmax  (t) 


A convenient  way  to  avoid  the  uncertainty  of  how  to  choose 
AN/&h  at  hmax  for  each  latitude  is  to  assume  that  the 
normalized  quantity  (aN/A^  ) ^/Nmax  ( t ) (a  profile  sha^e 
parameter)  does  not  change  appreciably  over  the  latitude 
range  covered  by  the  ray  path  in  the  F-region,  The  quantity 
may  therefore  be  evaluated  above  St.  John's  where  actual 
Nmax(t)  and  Nmax(t+/^t)  observations  are  available.  Thus: 


K*UN/&K)i4m.t  J 
Nma M 


ttiAtldt 

Nmax  It) 


wlhmax)  Jstau'i 


?c  w 


and  equation  (4.16)  becomes,  for  all  other  latitudes: 

Nmax(t*At)  ■ Nmax(t)(l  ♦ Cw(hmax))  (4.17) 
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We  now  have  only  one  adjustable  parameter  for  the 
distorted  array,  the  constant  K in  (4.14).  This  constant  is 
actually  a function  of  magnetic  and  electric  fields,  and  the 
time  (from  4.13).  We  choose  K in  order  to  reproduce  the 
observed  flat  TEC- i ncr eas i ng  Nmax  effect.  After  increasing 
Nmax  to  the  desired  value,  we  only  need  to  decrease  N(70U) 
and  N(1000)  by  sufficient  amounts  to  keep  the  Faraday 
rotation  constant.  For  the  distorted  array,  yzero  was  found 
from  available  ionosonde  data,  while  hmax  was  obtained  by 
decreasing  the  hinax  values  in  the  original  array  by  an 
average  of  20km  with  the  size  of  the  decreases  proportional 
to  w.  This  was  done  because  it  was  observed  that  hmax  tended 
to  drop  during  the  effect. 

On  occasion,  i.e.,  the  extreme  northern  and  southern 
ends  of  the  array,  use  of  equations  (4.14)  resulted  in 
unreasonable  values  for  N(700)  and  N(IOUO).  For  example, 
negative  values  for  electron  density,  or  N ( 700 ) <N( 1000 ) 
might  occur.  At  whatever  latitudes  this  happened,  the  values 
at  N(1000)  and  N { 7 0 0 ) there  were  obtained  by  distorting  the 
original  array  by  the  same  percentage  as  was  done  for  the 
nearest  latitude  point  at  which  reasonable  values  were 
obtained.  Unreasonably  high  values  of  Nmax  also  were 
obtained  using  (4.14)  at  the  extreme  poleward  end  of  the 
array.  These  points  are  north  of  the  ray  path  (except  when 
shifted  far  south),  so  they  will  not  contribute  to  Q.  In 

northernmost  reasonable  value  was  used 
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for  these  points  also. 

Having  discussed  the  approach  used,  let  us  now  present 
the  results  of  the  computer  model  simulations  of  ionospheric 
vertical  distortions.  In  Table  IV-XV  the  reader  will  find 
Nmax,  TEC,  and  Faraday  rotation  computed  for  both  Hamilton 
and  Goose  day  for  the  events  of  both  December  9 and  December 
15,  1971.  The  value  of  K used  in  equation  (4.14)  for  both 
dates  is  also  shown. 

Table  IV-XV 

Vertical  Distortion  Results 


STATION  Nmax(St . John's) 

TEC 

o 

K 

( 10*e  1 / on* ) 

(10u 

(deg) (meters ) 

el  /cina ) 

Date  Dec.  9, 

1971 

Or i gi nal 

Haml 1 ton 

3.95 

157 

Array-4:20  UT 

Goose  Bay 

5.4 

2.43 

108 

Di stor ted 

Hami 1 ton 

3,30 

131 

s 

Array-6 :00  UT 

Goose  Bay 

7.3 

2.4  0 

107  1 

.2x10 

Date  Dec.  15 

, 1971 

Original 

Hami 1 ton 

5.54 

219 

Array-2:25  UT 

Goose  Bay 

8.7 

3.04 

135 

DI stor ted 

Hami 1 ton 

6.81 

270 

c 

Array-b:45  UT 

Goose  Bay 

18.7 

3.04 

133  1 

. 12  x 10s 

The  values  of  K we  have  used  to  reproduce 

the  observed 

flat  TEC- i ncr eas 1 ng  Nmax 

effect  can  be 

used 

to  find  an 

approximate  value  for  the 

dawn-dusk  electric 

field. 

From 

(4.13)  we  arri ve  at: 

Em' 

-JLEks 

COS(ftLT)*t 

tf.lS) 

where  LT  Is 

In  hours,  E# 

o is  in  volts/meter. 

B^s 

is  in 

webers/meter*-= 

(vo  1 ts ) (sec)/metera’,  K Is  In  meters,  and 

At  1 s 
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in  seconds. 

In  (4.18)  we  use  the  times  midway  between  the  times  of 
the  original  array  and  of  the  disturbed  array  shown  in  Table 
IV-XV  for  the  local  time,  the  time  intervals  shown  in  Table 
IV-XV  for  At#  arid  the  value  of  .303x10**  wb/inx  for  Be^.  The 
calculation  gives  E*>*- . 67nV/m  for  the  Dec.  9 event,  and 
E*»=-.23mV/m  for  Dec.  15.  We  will  discuss  these  numerical 
E»*  values  further  In  sec.  4.4  when  we  take  up  the  question 
of  the  mechanism  for  the  flat  TEC- i ncreas i ng  Nmax  effects. 

The  distorted  arrays  for  Dec.  9 and  Dec.  15  are  graphed 
in  figures  I V—  2 5 and  I V- 24.  Direct  comparison  with  figures 
IV- 16  and  IV-16  reveals  the  type  of  distortions  that  have 
been  invoked.  The  dramatic  decreases  in  N(100d)  and  N(700) 
are  clearly  evident,  as  are  the  increased  Nmax  values. 

The  distortion  on  Dec.  15  seems  to  be  more  severe,  in 
that  N(700)  and  N(10J0)  are  reduced  to  values  typically  less 
than  10*el/cm*  for  high  latitudes  where  Nmax  is  multiplied 
by  a factor  of  from  3 to  5.  The  reason  for  the  larger 
distortion  on  the  15th  is  clear  from  Table  I V— VIII.  The 
large  percentage  increase  in  Nmax  (97%)  necessitated  the 
greater  compensati ng  decreases  in  N(700)  and  N(100Q). 

Notice  that  even  with  this  large  distortion,  the  17% 
decrease  in  TEC  has  not  been  reproduced.  From  Table  IV-XV  we 
see  that  this  large  distortion  was  accompanied  by  constant 
TEC  (3.04xl0,a'e1/cm*’).  To  produce  a decrease  in  TEC,  even 
greater  distortions  would  be  required.  For  December  9, 
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however,  the  observed  changes  have  been  reproduced  in 
detail.  The  relatively  small  increase  in  Nmax  evident  in 
Table  I V—  V I I (31*)  can  be  accounted  for  by  a smaller 


vertical  distortion  than  in  the  December  15  case. 


One  effective  method  for  demons tra t i ng  how  these 
vertical  distortions  are  capable  of  producing  a flat 
TEC-increasing  Nmax  effect  is  to  graph  Ne(h)  along  the  ray 
path  up  to  1000km  both  before  and  after  the  Nmax  increase. 
This  is  done  for  Dec.  9 in  figure  IV-25,  and  for  Dec.  15  in 
figure  1V-26.  The  anomalous  increase  in  the  distorted 
array's  value  of  Ne  at  about  170km  height  in  figure  IV-25  is 
due  to  the  lower  value  of  yzero  a^  the  highest  array 
1 at  i tu'-'es . 

One  can  see  from  these  figures  that  the  large  increase 
in  electron  concentrations  around  the  F-2  peak  and  including 
the  bottomside  is  compensated  for  by  lower  values  of  Ne  for 
the  distorted  array  from  about  40Ukm  out  to  the  satellite. 

As  previously  noted  (sec.  2. 1.3, 3),  the  computer  model 
adds  up  the  vertical  TEC  above  each  latitude  in  the 
interpolated  array.  This  vertical  TEC  is  directly  added  up 
to  1000km,  and  then  extrapolated  to  20u0km  using  a scale 


height  of  1000km,  as  explained  in  section  2.2.1.  The 
vertical  TEC  up  to  2 )00km  can  then  be  compared  with  the 


In  Table  IV-XVI, 


changes  in  the  vertical  TEC  up  to 


Ne(h)  Along  Ray  Path 
DEC.  9,  1971 


Original  Array 

Distorted  Array 


Ne(h)  Along  Ray  Path 

DEC.  15,  1971 

Original  Array 

Distorted  Array 


4 6 8 10  12 

Electron  Density 
(104el/cm3) 
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2000km  above  St.  John's  are  compared  with  corresponding 
changes  In  the  equivalent  vertical  TEC  at  Goose  Bay  for  the 
distorted  arrays  used  to  simulate  the  Dec.  9 and  Dec.  IS 
events.  Notice  that  the  values  of  equivalent  vertical  TEC 
are  all  higher  than  the  corresponding  vertical  TEC  values. 
This  confirms  our  conclusion  of  sec.  2.2.4  that  under  winter 
nighttime  conditions  Faraday  rotation  measures  the  total 
electron  content  of  a vertical  column  up  to  a height  greater 
than  2000km. 

It  is  encouraging  to  notice  that  the  changes  In 
vertical  TEC  are  only  about  12%/  because  this  does  not 
contradict  seriously  the  hypothesis  that  the  Faraday 
rotation  measured  at  Goose  Bay  monitors  changes  in  the  TEC 
of  a vertical  column  near  St.  John's. 


Table  IV-XVI 

Vertical  vs.  Equivalent  Vertical  TEC  - Distortion  Method 

(101*-  el /cm*) 


Dec.  9/  1971 


Original  Array 
Distorted  Array 
Dec.  15/  1971 


Original  Array 
Distorted  Array 


Vertical  TEC 
up  to  2 0 .) 0 km 
above  St.  John's 

2.10 

1.86 

Vertical  TEC 
up  to  2000km 
above  St.  John's 

2.  54 

2,86 


Equ i va 1 ent 
Vertical  TEC 
along  ray  path 
to  Goose  Bay 
2.43 
2.40 

Equi valent 
Vertical  TEC 
along  ray  path 
to  Goose  Bay 
3.04 
3.04 


Comparison  of  Table  IV-XV  with 
reveals  that  an  increase  in  TEC 


figures  IV-3  and  I V—  6 
at  Hamilton  has  been 
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reproduced  for  Dec.  15.  This  effect  increased  electron 
densities  at  lower  latitudes.  The  F-2  peak  height  is  at 
approximately  39  degrees  geographic  latitude  along  the  ray 
path  from  Hamilton,  a latitude  apparently  Inside  the 
plasmasphere  for  both  Dec.  9th  and  15th.  One  might  Initially 
propose  that  substantial  dumping  of  plasma  from  the  magnetic 
tubes  of  force  occurred  on  the  9th  but  not  on  the  15th  as 
this  would  explain  the  TEC  increase  observed  from  Hamilton 
on  the  9th,  which  we  were  not  able  to  reproduce  using  our 
vertical  distortion  technique. 

4. 3. 2. 4 Shift  and  Distort  Results 

As  pointed  out  by  Mendillo  (1973),  theory  predicts 
vertical,  longitudinal,  and  latitudinal  drifts  due  to  a 
dawn-dusk  electric  field.  While  we  have  not  been  able  to 
investigate  east-west  drifts,  we  have  looked  at  both 
vertical  and  latitudinal  motions.  One  might  ask  whether  a 
combination  of  vertical  and  latitudinal  drifts  could  account 
for  the  observed  flat  TEC- i ncreasing  Nmax  effect. 

To  investigate  this  question  we  have  utilized  the 
latitude  shift  capability  of  the  computer  model  in 

combination  with  the  distorted  arrays  plotted  In  figures 
IV-23  and  I V-24.  Results  are  shown  in  Table  I V—  XVII  for  0 
degree  to  -3  degree  shifts  on  Dec.  9 and  0 to  -u  degree 
shifts  on  Dec.  15.  One  should  compare  the  results  for  the 
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Table  I V-XV 1 1 
SHIFT  AND  DISTORT  RESULTS 

Date  Description  Station  O Model  Nmax 

(degrees)  4 2 0 -km  point 


I 

Dec.  9 Original  Array 

( 1J4  e 1 / cm5  ) 

Hami 1 ton 

157 

7.2 

1 1 

Goose  Bay 

103 

5.4 

1 

* 

Distorted  Array 

Hami 1 ton 

131 

B. 9 

1 

Shift  = 0 

Goose  Bay 

107 

7.3 

Distorted  Array 

Hami 1 ton 

107 

«.7 

Shift  *-l 

Goose  Bay 

155 

7.1 

K i 

i * 

Distorted  Array 

Hami 1 ton 

100 

OO 

• 

f i 

h 

1 a 

Shift  = - 2 

Goose  Bay 

247 

7.3 

1 1 

Distorted  Array 

Hami 1 ton 

102 

o . 6 

j 

Shift  —3 

Goose  Bay 

322 

8.8 

t 

J 


■L.. 


m 


u>  I WW: vm.  I j ipggppmui 


Date 


Dec. 
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Table  I V—  XVII  (continued) 

Description  Station  fl  Model  Nmax 

(degrees)  420-knn  point 
( 104  e 1 /cm4  ) 


Or i gi nal  - 

Ar  ray 

Hami 1 ton 

219 

15.3 

Goose  Bay 

135 

8.7 

D1  s tor  ted 

Ar  ray 

Ham i 1 ton 

270 

28.3 

Shift  = 

0 

Goose  Bay 

135 

13.7 

ni stor ted 

Ar  ray 

Ham i 1 ton 

264 

27.2 

Shift  = 

-1 

Goose  Bay 

167 

17.4 

Di  s to  r ted 

Ar  ray 

Hami 1 ton 

2 4 8 

26.2 

Shift  = 

-2 

Goose  Bay 

151 

16.2 

Ui stor  ted 

Ar  ray 

Hami 1 ton 

252 

25.4 

Shift  = 

-3 

Goose  Bay 

119 

la . 5 

Di stor ted 

Array 

Hami  1 ton 

213 

24.6 

Shift  = 

-4 

Goose  Bay 

194 

22.8 

Di  stor ted 

Array 

Hami 1 ton 

207 

2 3.9 

Shift  = 

-5 

Goose  Bay 

2 7y 

24.0 

Di  s tor  ted 

Ar  ray 

Ham i.l  ton 

186 

19.  5 

Shift  - 

• U 

Goose  Bay 

323 

38.  5 
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original  array  with  the  distorted  array  from  Goose  Bay  for 
each  shift. 

For  December  9 southward  shifting  greatly  increases  the 
Faraday  rotation  at  Goose  Bay,  and  is  thus  inconsistent  with 
the  observed  effect.  On  Deceiber  15,  the  -3  degree  shift 
gave  a 12 % decrease  in  Faraday  rotation.  This  compares  well 
with  the  observed  11%  decrease  (see  Table  I V—  VIII  anil  figure 
IV-6),  Remember  that  with  the  vertical  distortion  method 
alone,  a decrease  in  O could  not  be  realistically  achieved. 
These  results  imply  that  a combination  of  latitude  shift  and 
vertical  distortion  were  at  work  on  Dec.  15  but  only 
vertical  distortion  on  December  9. 

One  might  try  to  contend  that  to  reach  a -3  degree 
shift,  one  first  has  to  pass  through  the  -1  and  -2  degree 
shifts,  where  increases  in  O occurred.  However,  If  the 
original  array  represents  the  beginning  of  the  event,  and 
the  distorted  array  shifted  3 degrees  southward  represents 
the  peak  of  the  effect,  the  intermediate  stages  are 
represented  by  smal ler  distortions  than  that  of  the 
distorted  array.  In  other  words,  it  is  important  to  realize 
that  the  shift  results  * presented  in  Table  I V— XVI I do  not 
represent  a temporal  sequence,  but  only  a number  of  possible 
final  stages  of  the  two  events. 

intermediate  stages  with  smaller  distortions  just 
mentioned  have  not  been  Investigated,  as  this  question  opens 
up  a who  1 e new  level  of  complexity  in  the  calculations.  A 
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computer  program  which  automatically  distorts  input  arrays 
as  a function  of  time  could  be  developed  to  investigate  this 
question  further. 

O and  Umax  obtained  from  the  computer  model  for 
Hamilton  and  Wallop's  Island,  respectively,  are  also  given 
in  Table  I V- XVII.  These  results  should  be  compared  with 
figures  I V—  3 and  I V—  6 . Notice  that  the  increase  in  Cl  at 
Hamilton  has  been  reproduced  on  Dec.  15  using  shifts  of  -1, 
-2,  or  -3  degrees.  This  agrees  with  the  conclusion  that  a 
comoination  of  vertical  distortion  and  latitudinal  motions 
occurred  on  Dec.  15.  An  increase  in  O at  Hamilton  is  not 
obtained,  however,  on  Dec.  9.  This  is  not  a cause  for 
concern,  however,  since  we  have  already  concluded  that  a 
southward  shift  of  the  electron  density  vs.  latitude  array 
did  not  occur  on  Dec.  9. 

We  might  conclude  this  section  by  remarking  that 
although  model  results  seem  to  Indicate  that  a combination 
of  latitude  shift  and  vertical  drift  occurred  on  Oecenber 
15,  further  studies  are  needed  to  model  the  temporal 
evolution  of  the  event  on  that  night. 

4.3.3  Conclusions 


Computer  model  studies  were  undertaken  in  an  attempt  to 
simulate  the  flat  TEC- I ncreas i ng  Nmax  effect  which  occurred 
during  the  nights  of  Dec.  8-9  and  Dec.  14-15,  1971.  These 


studies  indicate  that  an  equatorward  shift  of  the  Ne(h) 


pattern  existing  before  the  effect  cannot  by  itself  account 
for  the  observations  on  either  night.  A vertical  distortion. 


responsible  for  the  effect  on  Dec.  9,  and  can  produce  a flat 


TEC  combined  with  an  Increase  in  Nmax  on  Dec.  15.  A 


combination  of  latitude  shift  and  vertical  distortion  can 


reproduce  the  observed  decrease  in  TEC  on  Dec.  15,  but 


cannot  reproduce  the  obser vat  I ons  on  Dec.  9 
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4,4  Mechanism  for  the  Slab  Thickness  Decreases 


It  Is  clear  that  the  vertical  distortions  of  the  winter 
nighttime  ionosphere  we  have  observed  near  L=4  are  substorm 
related.  On  all  five  nights  depicted  in  figures  IV-1  to 
IV-6,  substorms  were  observed  at  Mar ssarssuaq,  the  nearest 
of  the  eleven  geomagnetic  observatories  whose  records  are 
combined  in  the  AE  index.  For  the  case  of  Dec.  9 (figure 
IV-3),  the  beginning  of  the  distortion  of  the  ionosphere 
above  St.  John's  is  coincident  with  the  substorm  which 
occurred  at  4 UT  at  Mar ssarssuaq.  In  the  other  four  cases  a 
substorm  was  observed  at  Narssarssuaq  earlier  in  the  course 
of  the  night,  while  the  ionosphere  above  St.  John’s  was 
still  decaying.  There  could  have  been  an  immediate  response 
to  the  substorm  which  was  masked  by  the  underlying  decay. 
This  possibility  is  sketched  schematically  in  figure  I V—  2 7 . 


t\lmax 


figure  IV- 2 7 
The  areas  shaded 


are  the  same. 


time 


In  figure  IV-27,  possible  responses  of  the  ionosphere 
to  substorms  occurring  earlier  during  the  night  (while  the 
ionosphere  Is  still  decaying)  and  near  midnight  (when  the 


l ....  
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ionospheric  densities  have  reached  a state  of  stable 
equilibrium)  are  depicted.  Clearly , the  Increased 
concentration  near  midnight  would  appear  as  a "bump"  while 
that  occurring  earlier  would  not.  From  this  kind  of  argument 
it  becomes  apparent  that  immediate  response  by  the 
ionosphere  near  St.  John's  to  the  substorms  observed  at 
Marssarssuaq  could  have  occurred,  even  though  the  appearance 
of  such  a response  would  not  be  obvious  in  the  data. 

It  is  true,  however,  that  on  four  of  the  five  nights 
the  greatest  part  of  the  increases  in  St.  John's  peak 
density  occurred  a few  hours  after  the  substorms  occurred  at 
Nar ssar ssuaq,  at  times  when  the  ionosphere  would  be  expected 
to  reach  a stable  level  of  Ionization  during  magnetically 
quiet  conditions.  The  magnetic  conditions  were  not  quiet. 


however. 

on  those 

nights,  so 

one 

might  speculate 

that 

substorms 

might  have 

occurred 

near 

St.  John's  and 

near 

mldni ght 

whi ch  wer  e 

not  picked 

up 

at  Narssarssuaq. 

Even 

small  disturbances  might  greatly  affect  the  low  levels  of 
ionization  present  at  those  times. 

Akasofu  et.  al.  (l‘J/3)  showed  that  the  stations  whose 
magnetic  records  comprise  the  AE  index  would  not  record 
substorms  occurring  north  of  the  auroral  zone  along  a 
contracted  auroral  oval,  l.e,,  during  quiet  times.  Thus,  if 
a substorm  occurred  during  relatively  quiet  times,  its 
associated  enhanced  dawn  to  dusk  electric  field  might  cause 
ionospheric  distortions  without  a magnetic  disturbance  being 
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recorded  at  Narssarssuaq,  or  other  AE  stations.  One  may  also 
argue  that  the  11  magnetic  stations  whose  records  comprise 
the  AE  index  do  not  adequately  cover  all  longltiudes,  so 
that  a substorm  occurring  near  the  longitude  of  St.  John's 
might  not  be  detected  at  Narssarssuaq,  8.6  degrees  of 
magnetic  longitude  to  the  east.  Allen  ( ly 75)  has  confirmed 
that  this  might  be  a possibility,  although  Allen  and  Kroehl 
(1973)  concluded  that  the  11  stations  provide  adequate 
longitudinal  coverage  to  pick  up  most  substorms.  In  short, 
there  is  certainly  no  evidence  for  the  contention  that  all 
substorms  which  affect  the  ionosphere  near  L-4  are  reflected 
in  the  AE  Index,  We  therefore  can  contend  that  the 

distortions  which  occurred  on  all  five  nights  were 
associated  with  magnetic  disturbances. 

The  mechanism  by  which  a substorm  affects  the 

Ionosphere  near  L=4  is  the  enhanced  convection  electric 
field  associated  with  it.  This  enhanced  dawn-dusk  electric 
field  has  radial  and  tangential  components  which  map  along 
magnetic  lines  of  force  from  the  equatorial  plane  to  the 
ionosphere.  The  radial  component  causes  east-west  drifts, 
while  the  tangential  component  produces  both  vertical  and 
north-south  drifts.  While  we  have  not  been  able  to 
investigate  possible  effects  of  east-west  drifts  on 

ionospheric  density  profiles,  we  have  shown  that  both 

vertical  and  north-south  drifts  are  important  in  producing 
some  of  the  observed  effects. 
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There  are  fields  other  than  that  due  to  convection 
which  do  contribute  to  the  total  electric  field  in  the 
ionosphere.  These  are  the  corotation  electric  field  and  the 
dynamo  electric  field.  In  a review  paper  dealing  with  these  i 

various  electric  fields,  Mozer  (1973)  concluded  that,  on  the 
average,  the  ionospheric  electric  field  near  L=4  is 
primarily  of  magnetospher ic  origin  whether  the  pi asmaspher e 
Is  equatorward  or  poleward  of  L*4.  Under  magnetically 
disturbed  conditions,  the  convection  electric  field  is 
greatly  enhanced.  Table  IV-XVIII,  reproduced  from 
(Carpenter,  1970),  shows  how  the  dawn-dusk  electric  field  is 
enhanced  during  a moderate  magnetic  storm  period. 


Table  IV-XVI I I 
from  (Carpenter,  1970) 

Moderate  Magnetic  Storm  Period  Conditions  of 

prolonged  quiet 

Peak  Values  Typical  values  Average 

during  substorms  between  level  mV/m 

substorms 

mV/m  mV/m  mV/m 

Dusk*  1-4  ~0.4-0. 6 — 1 -0.1 

Midnight?  0.5-2  0.05-0.1  -0.2-0. 3 <0.1 

*based  on  statistics  of  the  position  of  the  plasma  bulge. 
Tbased  on  tracking  of  whistler  paths. 
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From  the  computations  In  the  previous  chapter  we 
derived  values  of  0,t>7mV/m  for  the  December  9th  event  and 
0.23mV/m  for  the  December  15th  event.  The  0.67  value  is  in 
agreement  with  the  peak  values  at  midnight  during  substorms 
in  Table  I V— X VIII,  and  of  course  a substorm  was  observed  at 
the  commencement  of  the  increase  in  Nmax  on  December  9.  The 
0.23mV/m  result  for  December  15  is  greater  than  typical 
values  between  substorms  in  the  table#  and  fits  well  into 
the  average  level  for  a magnetic  storm  period  near  midnight 
in  that  table.  This  seems  to  support  our  contention  that  the 
Nmax  increases  we  have  observed  are  due  to  substorm  related 
enhancements  of  the  convection  electric  field. 


i 

j 

j 
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One  final  comment  should  be  made  concerning  figure 
IV-27,  which  has  already  shed  light  on  the  relation  between 
substorms  and  nighttime  Nmax  increases.  This  figure  also  is 
most  relevant  to  the  question  of  why  the  nighttime  Nmax 
Increase  effect  is  primarily  a winter  phenomenon.  Only 
during  the  winter  months  is  Nmax  relatively  flat  for  many 
hours  during  undisturbed  nighttime  conditions#  at  which 
times  a ‘'bump'*  in  Nmax  would  be  obvious  In  the  data.  For  the 


same  reason#  most  of 
February  and  March 
classified  as  "low 
instance#  large  bumps 


: * 


the  bumps  observed  during  October, 
necessarily  fall  into  the  group  we 
bumps,"  while  during  December,  for 
can  occur  with  much  greater  frequency. 
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